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PREFACE 


It is widely felt that when the equations of mathematical 
physics are written out in full Cartesian form the structural 
simplicity of the formulae is often hidden by the mechanical 
labour of writing out every term explicitly. Attempts have 
been made to reduce this labour by one form or another of 
vector algebra; but it has always seemed to me that this 
method both introduces new difficulties and is insufficiently 
general, Thus the product of two vectors, in vector lan- 
guage, means one of two things, either the scalar or the 
vector product, and it is not physically obvious why just 
these functions of the vectors should arise and no others. 

The use of tensor notation, with the summation con- 
vention, carries out as great a simplification of the writing 
as does vector notation. The notation has actually at- 
tracted attention owing to its applications in the theory of 
relativity, but for ordinary purposes two great abbrevia- 
tions may be made. We use rectangular Cartesian axes; 
the result is that the distinction between covariant and 
contravariant vectors disappears, and with it the terms 
arising from curvature of the surfaces of reference. The 
formidable character of most of the formulae of the theory 
of relativity is absent from the formulae of tensors referred 
to Cartesian axes, The tensor method is a necessity for re- 
lativity; for applications in dynamics, electricity, elas- 
ticity, and hydrodynamics it is a great convenience. 

It is found that the scalar and vector products are not 
the only functions of two vectors that arise, though the 
theory provides reasons why they are important in many 
applications. There is also a symmetrical product, which 
ordinary vector notation is completely unable to express. 
In tensor notation it arises naturally as a symmetrical 
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tensor of the second order. The system of moments and 
products of inertia of a rigid body constitutes such a 
tensor; so do the stress components and the strain com- 
ponents in an elastic solid. 

The present method, like vector notation, is of use prin- 
cipally in proving general theorems. In concrete applica- 
tions there is usually some asymmetry about the coor- 
dinates that makes it necessary to abandon the tensor 
form at some stage in the work. It has been said that 
vector equations are like a pocket map, and it has been re- 
plied that a pocket map has to be taken out of the pocket 
and unfolded before it is of any use. The same applies to the 
tensor method, and for the same reason; but it has the 
great advantage that it is not a new notation, but a concise 
way of writing the ordinary notation, so that the unfolding 
can be carried out more conveniently when occasion arises. 

What is usually called Statics is treated in Chapter v, 
after Dynamics. I consider this to be the proper order, be- 
cause Statics is a special case of Dynamics, and many of its 
formulae have physical significance for reasons explained 
in Dynamics. The customary reversal of the order is due, 
I believe, to the fact that an introduction to mechanics has 
to be given at schools before the students have received 
any training in calculus; but this need not influence 
students working for a university examination. 

It should perhaps be stated that the object of this work 
is to illustrate the use of tensor methods; it does not claim 
to give a complete theory of all the subjects touched, re- 
ference for which must be made to the standard text-books. 

I must express my gratitude to Mr M. H, A. Newman, 
Miss L, M. Swain, Dr 8. Goldstein, and Dr Bertha Swirles 
for assistance at various stages in the work, and to the 
staff of the University Press for their care in the printing. 


HAROLD JEFFREYS 
September 1931 


NOTE 


Since this book was written most of the material in it 
has been incorporated in Chapters 2 and 3 of Methods of 
Mathematical Physics, by my wife, Bertha Swirles Jeffreys, 
and myself. The present reprint has been made partly 
because of a continuing demand for a treatment of Car- 
tesian tensors by themselves; partly because some results, 
notably on the thermodynamics of an elastic solid and the 
circulation in viscous flow, are not given in textbooks of 


the special subjects. 
HAROLD JEFFREYS 
January, 1952 


NOTE ON THE SEVENTH IMPRESSION 


In the present reprint a few small changes have been made. 
In Ex. 3, p. 15, the components are given explicitly ; these 
have been found useful in developing the strain energy for 
an elastic sphere under rotation. Omissions have been 
corrected on pp. 81 and 82. As these do not affect the 
later work they appear to have given readers no trouble. 


HAROLD JEFFREYS 
December, 1968 
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CHAPTER TI 
CARTESIAN TENSORS 


If we have two sets of rectangular axes (Ox, Oy, Oz), 
(Oz', Oy’, Oz') at the same origin, the coordinates of a 
point P with respect to the second set are given in terms 
of the coordinates with respect to the first set by the 
equations 

a! = Le + my + me 

y = he + my + m2 (1). 

a! = ha + mgy + Nye 


The quantities (1,,m,, 7, -..,% ) are the cosines of the 
angles between the various axes; thus /, is the cosine of the 
angle between the axes Ox’ and Oz, n, is the cosine of the 
angle between Oy’ and Oz, and so on. It follows that the 
coordinates (z, y, 2) can be expressed in terms of (x’, y’, 2’) 
by the relations 

w= he’ + hy’ +he 

y= me’ + myy! + Mg2’ (2). 

Z= Me + Nyy! + Ng" 


We can shorten the writing of (1) and (2) considerably 
by a change of notation. Instead of (a, y, z) let us write 
(1, %g, %g), and instead of (x, y’, z') write (2', 2’, 2%’). We 
can now say that the coordinates with respect to the first 
set of axes are w,, where i may be 1, 2, or 3; and those with 
respect to the second set are 2,’, where j may be 1, 2, or 3. 
Then in (1) each coordinate 2, is expressed as the sum of 
three terms depending on the three z,. Each 2, is associated 
with the cosine of the angle between the direction of that 
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x, increasing and that of x,’ increasing. Let us denote this 
cosine by a,;. Then we have, for all values of j, 


yf = yyy + MysXy + Aas 


- ney 37m (3). 


Conversely (2) can be written 
- , 

x patois 70% (4), 
the a,; having the same value as in (3), for the same values 
of ¢ andj, because it is in both cases the cosine of the angle 
between the directions of x, and 2,’ increasing. 

In mathematical physics we often have to deal with sets 
of three quantities in relation to a set of axes, of the general 
form u, (that is, ,, %~, Us), and such that in relation to a 
different set of axes the corresponding quantities are 
(tty', Ua’, Us’), which satisfy the relations 


uy Seen eaee (5) 
and u= ES ayn’ (6). 


j-1,2,3 
Such sets of three quantities are called tensors of the first 
order, or vectors. The individual w,, vu, v3 may be called the 
components of the tensors. 
Clearly if we multiply all of the u, and wu,’ by the same 
quantity m we get 
mu, = ae sat (mu) (7), 


8o that mu, is another tensor of the first order. 
Again, if we have two tensors of the first order, u, and 
v,, we shall have 
us +o! = , ea 119 (te + %) (8), 


so that u, + v, is a tensor of the first order. 
We notice that each of the equations (3) to (8) is really 
a set of three equations; where the suffix i or j appears on 
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the left it is to be given in turn all the values 1, 2, 3, and the 
resulting equation is asserted in each case. In each such 
equation the right side is the sum of three terms, obtained 
by giving) or i the values 1, 2, 3in turn and adding. Wher- 
ever such a summation occurs a suffix is repeated in the 
expression for the general term; where there is a summa- 
tion for all values of j the general term, such as a,,%,, con- 
tains j twice. We make it a regular convention that, unless 
the contrary is stated, whenever a suffix is repeated it is 
to be given all possible values and that the terms are to be 
added for all. Thus we write (5) as simply 


Uy = ys, (9), 


the summation sign being automatically understood by 
our convention. Then (9) really means three equations, 
with three terms on the right of each, but we can by means 
of our conventions express all of the twelve terms com- 
pactly by the single equation (9). 

There are single quantities, such as mass and distance, 
that are the same for all sets of axes. These are called 
tensors of zero order, or scalars. 

Consider now two tensors of the first order, u, and v,. 
(When we write “a tensor u,” we mean of course a tensor 
of the first order whose components are t, U,, U3. This is 
another piece of shorthand.) Suppose each component of 
the one multiplied by each component of the other; then 
we obtain a set of nine quantities expressed by u,v,, where 
each of i and & is independently given all the values 1, 2, 3. 
The components of u,, v, with respect to the other set of 
axes are u,’, v;' say; and 

Uy Dy! = (eye) (AerYe) 
= UA Uy (10). 


The suffixes i and & are repeated on the right. Thus (10) 
represents nine equations, each with nine terms on the 
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right. Each term on the right is the product of two factors, 
one of the form a,;a,:, depending only on the orientation of 
the axes, and the other of the form u,v,, representing the 
products of the components referred to the original axes. 
In this way the various u,’v;/ can be obtained in terms of 
the original u,v. But products of two vectors are far from 
being the only quantities satisfying this rule. In general a 
set of nine quantities wy, referred to a set of axes, and 
transformed to another set by the rule 


We! = Ag Ay Wen (11), 


is called a tensor of the second order. 

We may go on similarly to construct and define tensors 
of the third, fourth, and higher orders. Thus a set of quan- 
tities that transforms like 2,2,%,,2,... is called a tensor of 
order n, where n is the number of factors in this product. 

When we say that a certain set of quantities is a tensor 
of any order , we mean that we have ways of specifying 
its components with respect to any set of axes, and that the 
components with regard to any two different sets of axes 
are related according to the rule appropriate to tensors of 
that order, and in particular to the products of the coor- 
dinates with n factors. For instance, if we say that u, is a 
tensor of order 1, we are not simply defining w,’ as meaning 
aju,;. We are supposing both that u,’ has a meaning, such 
as a displacement or a velocity, with reference to the axes 
of z,', and that the value of each component is equal to 
a,;u,. Thus the statement that any set of quantities is a 
tensor is not a mere convention, but a statement capable 
of test and therefore needing proof. In (7) and (8), for 
example, our data are that wu, and u,’ are the components 
of a vector with regard to two different sets of axes. We 
prove that the sets of quantities obtained by multiplying 
both by the same quantity are related according to the 
vector rule; and therefore the products are vectors. 
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If we interchange j and 1 in (11), we get 
Wy! = Aarts We (12). 
But on the right i and k are “dummy suffixes”; that is, 
they are to be given all possible values and the results 


added. It is unimportant which of them we call i and which 
k; we may therefore interchange them and get 


Wy! = Aeris Wey = Ay Aner (18). 
‘Thus w,; transforms according to the same rule as wy, and 


therefore is another tensor of the second order. The im- 
portance of this is that if we know the set of quantities 


arranged 
Wy Wie Wis 
(« Weg vs) (14) 
Wy Wyq Wa, 
to be a tensor of the second order, then the arrangement 


Wy Uy Wy 
( Wye Wy Wye ) (15) 
Wyy Wey Way 
is another tensor of the second order. Therefore the sets 
(wa, ++ Wye) and (wy, — W,;) are tensors of the second order. 
The first of these has the property that it is unaltered by 
interchanging ¢ and &, and is therefore called a symmetrical 
tensor. The second has all its components reversed in sign 
when i and & are interchanged, and is called an antisym- 
metrical tensor. Clearly in an antisymmetrical tensor the 
“Jeading diagonal” components, i.e. those with i and k 

equal, are all zero. Also, since 


Wa =F (Win + Wes) +4 (Wa — Wer) — (18), 
we can consider any tensor of the second order as the sum 


of symmetrical and antisymmetrical parts. 
The gradient of a scalar is a vector. For if U is a scalar, 


6 CARTESIAN TENSORS 


its gradient is dU /dx, or 6U/ex, according to the set of 
axes. But 

OU _ ax BU OU 

Ox) Ox,’ Ox, Oa, 
so that the gradients transform according to the vector 
tule. Similarly the gradient of a vector is a tensor of 
order 2. For if u, and u,’ are the components of a vector 
with respect to two sets of axes, 

, , a 

ous saad s oe _ Ont oa (a5 U4) 


(17), 


Ou, 
= AsO aang (18), 


so that the rule of transformation is as in (11). 

Since 2, is a vector, it follows that @x,/02, is a tensor of 
the second order. But éz,/02, is unity if i = k and zero if 
i#k. Hence the set of quantities 5, such that 

bu= bn baa 1, 
Bie = Bip = Ban = 8x3 = 8 = On = 0, 
constitutes a tensor of the second order. We can prove this 
directly; for if we apply (11), 5,’ in the new system of 
coordinates should be given by 
85 = Ay tudin (19). 
The suffix & has to take all values 1, 2, 3. But if k#i, 
8, is 0, and the corresponding term is zero. If k=i, 
8 = 1, and the result of the summation with regard to k is 
By! = ayaa (20). 
But the a,, are the direction cosines of the axis of x,’ with 
regard to the a, and the ay, are those of x,’ with regard to 
z;. Hence a,;,1 is the cosine of the angle between x,’ and 


2;', and is equal to 1 if the axes are identical and to 0 if they 
are perpendicular. It follows that the result of the trans- 
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formation is that 5,,' = 1 if j = 1, and 8,’ = 0 ifj #1. It 
follows that the set of quantities 


10 0 
(° 1 *) (21) 
\ i ies 


is transformed into itself by the rule (11) and therefore is 
a tensor of the second order. 

If u, is a vector and we form the product 8,,.%,, we have 
a tensor of order 3. But now put m =k and add for all 
values of k. Since 8,=0 except for & = i, the only term 
different from zero is that for k = i, and this is w,. Hence 

Sate = Uy (22). 
This operation therefore replaces the suffix k by i. The 
tensor 5,, can therefore be called the substitution tensor. 

In the tensor wy let us put & =i, and in accordance 
with our convention add for all values of i. Then the 
corresponding quantity w,,’ is got by putting 1=j and 
summing; but 

Wy! = Aisa = 3a We 
= wy (23). 
Thus w,, transforms into itself and therefore is a scalar. 

‘This operation of putting two suffixes in a tensor equal 
and adding accordingly is known as contraction. In general 
it gives a new tensor, whose order is less by 2 than that of 
the original tensor. If for instance we contract the tensor 
Uv, we obtain 

UjVy = Uy Vy + UgYy + Wyre (24), 
which is the scalar product of u, and v,. 

Similarly the tensor u,v, yields the symmetrical and 
antisymmetrical tensors (u,v, + %,0;) and (uv,—u%v,). We 
may call these the symmetricaland antisymmetrical products 
of u, and »,. 
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The tensor du,/2, gives similarly, on contraction, a scalar 
Duy _ Orly , tly , Oly 
Be, Be, * By * Bay big 

which is known as the divergence of u,; while it gives also 
symmetrical and antisymmetrical tensors 
tly, Ou Ou, _ Ou, 
Ox, + i, end Oa, Oa,” 
The former has important applications, especially in the 
theory of elasticity and hydrodynamics; the latter is 
known as the curl or rotation of u,. The vanishing of the 
curl is the condition that w, may be the gradient of a scalar. 

All the above considerations can be extended to any 
number of dimensions. In n dimensions a tensor of order r 
has n* components. A tensor of order 2, in particular, has 
n® components. If it is antisymmetrical, the n diagonal 
components are zero, and the others are equal and opposite 
in pairs, Hence an antisymmetrical tensor of order 2 has 
4n (n — 1) independent components. Ifn=1, 2, 3, 4,... in 
turn, this number is 0, 1, 3, 6,.... It happens that in three 
dimensions the number of numerically independent com- 
ponents of an antisymmetrical tensor of the second order is 
equal to the number of components of a vector. Actually it 
can be proved that with any vector we can associate an anti- 
symmetrical tensor of the second order, and conversely. 
‘This is not true in any number of dimensions other than 3. 

Since the a,, are the direction cosines with respect to the x, 
of three perpendicular lines, they are connected by six re- 


lations Ay? + dy? + Gy)? = 1 
yy" + gq” + yg” = 1 (26), 

yg? + Ay” + yg? = 1 

Aya (hyy + Mantes + Aang = 0 
Ayytyy + Mpg + Ayan = O (27). 

Ay, Ayy + Aye + O51 % = O 
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We notice that the second and third of (27) both contain 
(a4, 4), 43,). We may therefore solve them for the ratios 
of these quantities. Thus 


CE: ae Sarr as 
gg Mag — Ugg Qpg Ugg yg — AygM%gq Tyga — Aggy Myo 
(28), 


say. Substituting in the first of (26) we get 
1 = k {ayy (Any — sage) + nr (455412 — ra M9) 
+ Ay (443432 — Mp3 %y»)} 


= a—hlay dy ay (29). 
Cm 
Also ies: 
He? { (dy gq — gg aa)? + (ggg — Mya %g2)* + (442 — M25 %2)"} = 1 
(30), 


But we have a general identity 
(a? + 6? + c?) (a’* + b’* + 02) — (aa’ + bb’ + cc’)® 
= (be' — cb’)? + (ca’ — ac’)? + (ab’ — ba’)? (31). 
Hence 
Ke [(ays? + ag? + yg") (Qia® + Oaa® + Aya") 
= (Gazi + A2a%5 + A52%3)*] = 1 (32). 
But on account of the second and third of (26) and the first 
of (27) the expression in brackets is unity, and therefore 
k=+1 (33). 
For any given transformation the determinant in (29) is 
therefore equal to + 1. Evidently its sign is reversed if we 
interchange any two of the suffixes j, for this interchanges 
two rows of the determinant; so that the sign is a matter of 
the numbering of the axes. If we start with a rigid frame 
attached to the axes a, and rotate it continuously till it 
is attached to the axes a,’, all the a,, vary continuously and 
therefore the determinant cannot change from + 1 to — 1 


2 yer 
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or from — 1 to + 1. If then x, goes to 2’, 2, to 24’, and 23 
to a,', the determinant is initially 


10 Oj=1 (34), 
010 
001 
and therefore 
k=-1 (35), 


and the determinant formed by the a,, is always + 1. 

If we have a frame of axes (x,2,%;) we can turn it by a 
continuous movement so as to bring x, along the old x, 
2, along the old a, and a, along the old x,. In this case we 
have A 

Dy! = Wy, Ty = Hy, Ty = (36), 
and 
Gy = 0, dy = 1, dy = 0; G2=0, Aye = 0, Oya = " (37). 
%3 = 1, Ay = 0, Oy = 0 

The determinant of the a,, is therefore 

0 1 Oj=1 (38), 

oot 

TniOy) O 


as before. Any rotation of the axes that does not alter the 

cyclic interchange of suffixes 1231231... therefore 

leaves the determinant equal to unity, and therefore so 

long as we always use right-handed or always left-handed 

axes the determinant of the a,, is + 1. 

With this restriction 
yy = Ugg gg — Ty%9g3 A = Tis To — cd (39), 

xy = Ay3 M2 — M3 %ra 

and therefore every direction cosine is equal to its first 

minor in the determinant. 
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These relations are of course identical with those ex- 
pressed in the usual notation of solid geometry by 
J = myn — MyNy; M = Nyly — Nyly; My = LyMg — lyymy 
(40). 
Now suppose that u, is a vector, and consider the set of 


quantities 
=; 0 & = (41). 
—-% 0 «) 


UW —-%m 0 
Apply (11) to this, taking 7 to be the number of the row and 
& that of the column. We see that w, enters as w,, and as 
— Wy. Its coefficient in w,;’ is therefore ay,a3,; — d3)4,,, 
and in all 


Wy! = (a5 gx — ys Ma) Uy + (gs y4 — 451) Ug 
+ (Qys4421 — 253) Ug (42). 
This is obviously zero if j = 1. If j #1 and if the other axis 
perpendicular to =,’ and z,' is z,’, and jlnjln is a cyclic 
order, the quantities in brackets are equal to (aj, dan, Gan)» 
This is true if 1 immediately succeeds j in the order. If 1 
precedes j by one place the signs are reversed. Hence if 
j= land/= 2, orj = 2and/= 3, orifj =3 and/=1, 


Wy! = AjyUy = Up! (43), 
and in the alternative case 
Wy! = — Uy! (44). 
‘Thus wie, 0 wl — my (45), 
Fen, 
u’ —%' 0 


and is of the same form as (41). Thus with any vector we 
can associate an antisymmetrical tensor of the second order, 
Conversely with any antisymmetrical tensor of the second 
order we can associate a vector. 
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We can proceed alternatively by considering the set of 
quantities €jm, defined by the condition that if any two 
of i, &, m are equal the corresponding component is 0; if 
i, k, m are all unequal and in cyclic order, the component 
is + 1; if the order is not cyclic, the component is — 1. 
Let us see whether this is a tensor of the third order. If so, 
we should have 

jn = U5 QeG mn Eve 
= Qj Mg1 Fyn + Ag Mgr%in + 54114an 
= yyy, Ay — Mg 5Mg1%in — yj %1%n (46). 
Now if, for instance, j = J, the right side is clearly zero and 
€sm' = 0. Ifj, 1, are all unequal, the expression is 


sy Ay Aas 
Gy Ag, Agr (47), 
Un Fn An 


which is equal to 1 if jin are in cyclic order and to — 1 if 
not. Hence the set of quantities <,., is transformed into 
itself by the rule for transforming tensors of order 3, and 
therefore constitutes a tensor of order 3. This is called the 
alternating tensor. 

Now consider the product €jm%), Where w, is a vector. 
This is a tensor of the fourth order. If we contract it by 
putting p = m and summing we get a second order tensor 
Wa = €iymtm If i= land k = 2, the only value of m that 
makes ¢,, different from zerois 3, and then éam= + 1. Hence 

Wy = ty (48). 
Ifi=2andk = 1, mis 8; but 213 is the reverse of eyclic 
order and €y; = — 1. Hence 
Wy = — Us (49). 
Similarly we find that the elements of w,, are 


0 Uy — 
(- Uy 0 *) (50), 
Uy — th 0 
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so that the antisymmetrical tensor associated with a vector 
can actually be obtained from it by multiplying by um 
and contracting. 

Again, suppose that we are given a tensor of the second 
order wy, and that we form a vector u,, by multiplying by 
€am and contracting twice. We have, if m = 3, 


Us = €uaMe = Gaia + San = Wig — War 

Thus Un = Cem ir (51). 
If we, is symmetrical, this evidently gives zero. If it is 
antisymmetrical the components of u,, are numerically 
twice those of wy. 

On account of the intimate relation between the vector 
and the antisymmetrical tensor we shall habitually denote 
the tensor w,, of (41) by wy, 80 that 

‘Uy = Ugg = Usa = 03 the = Ua» Uae = as tee (62) 

gy = — May Mag = — ths tha = — Ue 
It will always be seen at once whether the vector or the 
tensor is intended, since the former has one and the latter 
two suffixes. 

If we have any three vectors u,, v;,,, and consider the 
BCAAT € jum UiVe Wm, We see that 

perm Ui Vg Wy = Uy Ug Wy + UyV gy + Uy Vy Wy 
Ug dy Wy — UyVgWy — Uy Ug Wy 
=] % th th (53), 
My Uy Us 
W, Wy Ws 
so that we have a concise way of writing the determinant 
formed by the components of three vectors. If any two of 
the vectors are parallel this scalar vanishes. 

In associating a vector with an antisymmetrical tensor 
of order 2 a sign convention clearly arises. We make the 
positive signs in (41) lie one place to the right of the leading 
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diagonal. If then we have two vectors u, and v,, their anti- 
symmetrical product is 4,7, — %v;, and in the associated 
vector we give the positive sign to wg, when & follows ¢ in 
the cyclic order. Hence the components of this vector are 
taken to be 


(gy — ge, UP, — Uy Ugs th Ve — Uy). 
This vector is perpendicular to both the original vectors; for 
Ay (thgVy — Ua) + Uy (gt — Ua) + Uy (thxYe — tat) = 0, 
 (Ua% — Maa) + Ue (Warr — UY) + Us (tM — tet) = 0. 
We call it the vector product of u; and v,, and can save 

writing by denoting it by [u, Um. 

Similarly with the antisymmetrical tensor - a we 
associate a vector so as to leave the sign unaltered when 
k follows i in the cyclic order. Thus the components are 


(oO, HM, Be MH) Gy 
zy de,’ Oey my? Bx, : 
This is often called “curl w.” 

The Tensor €ix.€nys- Since this tensor is the product of 
two third order tensors, once summed, it is a tensor of 
the fourth order, i, &, m, p being arbitrarily assignable. 
Evidently if i= k or m= >p, the corresponding component 
is zero. 

If i =m, the contribution from any value of s is zero 
unless also k = p, and then 

€ixs = Emps = + 1, 
and the component is + 1. 

If i = p, then no value of s gives a contribution unless 
k =m. Then one of ¢,,, and €,,5, is + 1 and the other — 1, 
and the component is — 1. Hence the components of the 
tensor are as follows. 
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If i= m, k = p, the component is + 1, unless ¢ also = k, 
i= p,k =m, the component is — 1, unless é also = k, 
i=k or m= p, the component is 0. 
These results apply also to the tensor 
8imSin — Seo5ems 


and therefore 


€as€mos = Simdey — 81092 (55). 


EXAMPLES 
1. If u;, uy’, u,” are the components of a vector with regard to 
three sets of axes, prove that the values of 1,” are the same as 
would be obtained by transforming first from u, to u,’ and then 
from uj’ to u,". 
2. Prove that 34 = 35 8am = 0- 
8. Show that 
Kerem = ba18025em — BacBomBce + FaxBomex 
— BaxBoi5em + Sam Boi8cx — Sam Boi Ser 
4. Prove that €taSnta ~ Bim Com Em = 8» 
5. Prove that its€mpu = Gatnfamn = €tutpume 
6. Prove that if u,, vg, Wp, are vectors, 
[ew [ety PD mn = thon (04104) — Yon (44204 )o 
than [% Wm = oem tiem 


7. if A(u)=| 4 the ts |» 
Yn Ms Un 
Us, Use Usa 
prove that 


eee 3 (4) = Cima Yan® 
setts Min = “mn & (tes 
6A (u) = C492 €imn Yt Ym Yan 
8. Use Ex. 7 to prove the rule for the multiplication of deter- 
minants 
ma “mn & (tu) A (v) = 6A (u) A (v) = 6A (wr), 
where (Wr)ip = Uri 


CHAPTER I1 
GEOMETRICAL APPLICATIONS 


The displacement from any point to any other obviously 
constitutes a vector. The distance between the points is a 
scalar. If a,, y; are the coordinates of the points and r the 
distance between them, 
r= (yi — x)? (1), 
the square on the right indicating the scalar product of the 
vector into itself. Also the quantities (y, — 2,)/r constitute 
a vector. 
If we take a fixed point a, and consider points given by 
asa +hr (2), 
where r is a variable scalar and the J, are constants such 
that 
Ig=1 (3), 
(a, — 4)? = 1? (4), 
so that 7 is the distance of x, from «,. If we take another 
point y, such that 


w=aths (5), 
(ye — a)? = 8* (6), 
(ye — %)* = (8-7)? (7), 


and therefore the distances between «,, z;, and y, are such 
that the sum of two of them is equal to the third. Thus the 
points are on a straight line; and (2) gives the equations 
of the line in terms of the parameter r. The J, are the 
direction cosines of the line. 
If we take two lines through @, given by 
a= a+ lr (8), 
Y= a+ m8 (9), 
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the distance between «, and y; is given by 
(ye — %4)* = (8m, — rh,)? 
= st + 72 — Orslym, (10). 
But this quantity is also equal to s* + r* — 2rs cos 0, where 
0 is the angle between the lines. Hence the angle between 
two intersecting lines is given by 
cos 8 = l,m, (11). 
If two lines have the same direction cosines they are said 
to be parallel. If two lines do not intersect we can take a 
line through any point on one of them parallel to the other; 
then this line is inclined to the first at an angle given by (11). 
We can then use (11) to determine a unique quantity 
associated with any two lines, which we may call their 
inclination, whether they intersect or not. 
If we have a line given by (2) and y; is a point outside it, 
the line joining a, and y, subtends a right angle at «, if 
(ya — os)® = (ae, — %)* + (Ye — 2%)? 
= 78 + (yy — a — 1h)" 
= (ye— @)® = 2rly (ye — ay) + 2712 (12), 


and therefore 
r=] (y— &%) (13). 


This gives the projection of the displacement y,— a, on 
the line. The foot of the perpendicular is 


ay + Ly = a + Lily (Ye — O%)- 


Evidently 7 in (13) will be the same for all points y; such 
that ly, is constant. Hence 


liys=8 (14) 


represents a plane perpendicular to the line. 
If we take two intersecting lines given by (8) and (9), we 
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can find the equation of the plane containing them as 
follows. If this plane is 

m= PD (15), 
this equation must be satisfied by x, and y, for all values of 
rands. Hence 


Nyt =p (16), 
nl, =0 (17), 
nm, =0 (18), 
and from (15) and (16), 
nN, (%—- a) =O (19). 


Then (17), (18), (19) are three homogeneous equations in 
the m,;, and can be consistent only if 
€inm (24 — 4) Lp Mp = 0 (20). 
This is the equation of the required plane. Also the n, are 
proportional to 
Edam la Mm = (IyiMg — Lymg, tym — Lym, bmg —1ym,) (21). 
But ng=1 (22), 
(dys — Tyrmg)* + (lym, — bymg)® + (lym, — Lam,)* 
= (L8 + be? +h?) (my? + mg? + mg*) — (my + lymy + Lyng)? 
= 1-—cos*@ 
= sin? @ (23). 
Thus sin 0 y= + €amlmMm (24). 
The ambiguity in sign corresponds to a general one in 
specifying the parameter r of a point on a line. If r in (8) 
is taken negative, we get a point on the line on the opposite 
side of a, from those given by positive values of r. But if 
we reverse both r and the J, we still keep 1,4 = 1, and we 
still have the same point. We may take either direction 
along a line to be that of r increasing; if we reverse the 
direction the signs of all the 1; are reversed for the same 
point. 
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For any point on (8), 

€ sem ®lm = €cem (Ce + Let) Im = Esa Gebm (25). 
This is a constant vector for all values of 7, and may there- 
fore be considered as a property of the line. We denote it 
by 1. Then we have six properties of the line given by 
1,, l. These are coordinates of the line. They are connected 
by two relations, 


Ig=1 (26), 
Ue = Cem lidelm = 0 (27). 
The J, have a geometrical interpretation. Thus if we con- 
sider the plane 
%yly — thy = 1’ (28), 


this plane passes through the line. Also if 2, = a, 2% = a, 
(28) is satisfied for all values of x,, and therefore if we take 
a line through @, parallel to the x, axis, (28) represents the 
plane through (8) and this line. Two such planes determine 
the line, and therefore the J, and J,’ together determine the 
line. 

If we have two non-intersecting lines given by 


m= a, + rh, (29), 
Ye = Be + 3m, (30), 
the line 4 = Be + rl, (31) 


passes through , and is parallel to (29). The plane including 
(30) and (31) is, by (20), 

€sam (2 — Bi) eM = 0 (32). 
This therefore represents a plane through (30) parallel to 
(29). The plane through (29) parallel to (30) is 

sem (2% — @:) eMm = 0 (33). 
The distance between these planes is the projection of the 
line joining any two points on them upon a line perpen- 
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dicular to both. If a line perpendicular to both has direc- 
tion cosines n,, the shortest distance d between the lines 
is therefore given by 
dsin 0 = (8, — a) n, sin 8 
= + (Bi— &) Ciemle Mm 
= + {eumBiltn — €ameileMn} (34). 


But € tem Bi Mm = — Mx" 5 Comm ile = bn! (35), 
and therefore 
dsin @ = + (— hm, — ml’), 
so that, apart from the ambiguity in sign, 
dsin 0 = l,m, + ml,’ (36). 
Thus the shortest distance is directly expressible in terms 
of the coordinates of the two lines. 
Now consider two intersecting lines 
B= at rhs yg = ay + amy (37). 
The area of the triangle formed by «,, x, y; is Irs sin 0. 
The projections of these points on the plane 2, = 0 are 
(0, ag, ag), (0, 2, 3), (0; Yes Ys) and form a triangle whose 
area is 
3 =h}l q a 
Oo br lr 
0 m8 mss 
= ds (ym; — lym) 
= dren, sin 0 (38). 
Thus the projections of a triangle, and therefore of any 
plane area, on the coordinate planes are in the ratios of the 
direction cosines of the normal to the planes. A plane area 
can therefore be treated as a vector whose components are 
proportional to the direction cosines of the normal. 


1 a a 
1 &% 2% 
lt Y% 
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If we have a line given by 

= a + rl, (39), 
and 8; is a point not on the line, let us suppose the point 8; 
turned through an angle @ about the line. The foot of the 
normal from £, to the line is given by 
r=, (By — a) (40), 
and therefore the displacement from the foot of the normal 
to B, is equal to 
Be — fe + Lil (Be — &x)} = Ue? (By — 44) — Lily (Be — ax) 
= ly {he (Bi — &) — Le (Be — %)} 
(41). 
The magnitude of this displacement, p, is given by 
P= (Be — a)? — {h (Be — a}? 

= U2 (B; — @%)* — Lh, (Bi — @) (Be — &) 

= 3 fh (Br — &%) — bs (Be — a)? (42), 
the } being needed because in the double summation each 
pair of values of the suffixes would occur twice. 

‘The plane through £, and the line is 

Ny (% — @) = 0, 
subject to nl, = 0, 
m, (By — @) = 0, 
and is therefore 
€ sam (2 — &) Ly (Bm — &m) = 0 (43), 
while the n, are proportional to € xml; (Bm — &») and there- 
fore equal to + ésmls (Bm — Gm)/P- 

If now we turn f, through an angle 6, it receives a dis- 
placement p (1 — cos @) along the normal to (39), and a 
displacement p sin @ along the perpendicular to the plane 
(43). If it goes to y,;, we have therefore 
Ye — Be = — (1 — 008 8) Ly {le (Bs — 4) — 1 (Be — )} 

Sin 8 €gmle (Bm — Gm) (44). 
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If 6 is a small angle and we neglect 6°, the displacement is 
simply 
Ye — B= + 8 €emle (Bm — Em) (45), 
and in particular 
Y1— Br = {19 (Bs— 4) — 1,8 (Bz —a2)} (48), 
and from the additive form of this equation we see that the 
displacement is the sum of those given by separate small 
rotations 1,6 about axes through e, parallel to the co- 
ordinate axes. Conversely, displacements due to small 
rotations about axes through a point can be added vectori- 
ally as if all were applied to the system in its original 
position, and give the same total displacement as if they 
were compounded into a single rotation about an axis by 
the vector rule. We still have, however, to establish a sign 
convention. We decide that @ is to be taken positive if a 
turn about the axis of x, is from 2, towards z,. This would 
make 
v1 — Bi= — 9(Bs— 2); Y2— B2=9(Bi—%) (47) 
with 1, = 1, = 0, , = 1. Hence 


4 — Be= mle? (Bm — mn) (48). 
If we write 
L,0=%, 
Ye — Be= €aem (Bm — tm) (49). 


For instance, if «,, = 0, we have 
11 — B= 283 — MPs; Y2— Pe = M8, — 7,833 
Ys — By = MB, — 8 (50). 
For finite rotations we return to (44) or (41), keeping the 
positive sign in the second term. We may transfer the 
origin to @, to save writing. The coefficient of 8, in y, — B; 
is 
ba = (1 — 008 0) 1jly — 8in Oeiemlm (51), 
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for k # i; if k = i, the coefficient is — (1 — cos 8) (1 — J;1,). 
Thus 
ba = (1 — 008 8) (ih, — 84) — Sin Pom (52), 


where cy is the antisymmetrical tensor corresponding to 


1,, namely ‘i 
0 = 
(- Led i) (53). 
4h -h 9, 


Thus the displacement is represented in general by by.B,, 
where b, is a tensor of the second order, expressed as the 
sum of symmetrical and antisymmetrical tensors. The anti- 
symmetrical part is seen to be of the first order of magnitude 
in 6 and the symmetrical part of the second order, 


EXAMPLES 
1. Given that the general quadric surface is 
S=}4g27,a, + Bx, + 0=0, 

prove that the locus of the mid-points of parallel chords is a plane, 
and find the condition for this plane to be perpendicular to the 
chords. 

2. Find the condition that the line 

lat p=0 

may touch the quadric 8. 


OHAPTER IIl 
PARTICLE DYNAMICS 


The essence of particle dynamics is that the second de- 
rivatives of the coordinates of a particle with regard to the 
time are equal to functions of the position and velocity of 
the particle with reference to neighbouring particles. The 
relations therefore provide a set of differential equations 
to determine the coordinates. The equations for any 
particle can be put in the form 
mii, = ZX, (1), 

where m is the mass of the particle, X, are the forces due 
to the other particles, and the summation is for all the 
other particles*. It is a matter of experiment that this 
form holds when the axes are a certain type of Cartesian 
axes, which we call non-rotating and unaccelerated, or in 
brief dynamical. 

If we take a different set of non-rotating axes at the 
same origin we have 

By = Ay, (2), 
and since the axes are not rotating the a,, are constants. 
Hence by differentiation 

By = agit, (3), 

Hy = ays, (4), 
and therefore velocity and acceleration are vectors. 

The force X, on the particle due to some other particle 
is measured by the contribution to mi, due to the other 
particle; that is, the part of md, that would disappear if 
the other particle was removed to an infinite distance. It 

* Cf, Jeffreys, Scientific Inference, Chapter vimt, for a fuller analysis. 
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is actually measured by the acceleration. If we have the 
three acceleration components of the particle 1 due to the 
particle 2, we can find the acceleration in the direction of 
x,’ by the formula (4), and the form (1) will become 


mij = DX; (5); 
provided that we define X,’ by the rule 
Xf = aX, (6). 


The meaning of the force in any direction requires de- 
finition in any case: if we define it by its relation to the 
acceleration component in that direction, it follows auto- 
matically by (6) that force is a vector. It follows that the 
sum of any number of forces, obtained by adding the 
components separately, is also a vector. But the practical 
importance of the idea of force arises equally from the fact 
that in many cases the force components are known once 
for all from experience as functions of the coordinates and 
velocities. The form (6) then says, as a general principle, 
that the forces in any direction are additive. 

If we form the contracted or scalar product of md, and 
=X, by the vector #,, we get 


mi, , = DX ty (7), 
the left side of which is 4 ma), (The square of course 


implies the product of #, by #, and therefore the summation 
for i = 1, 2, 3.) We write 

T= Ym = Ym (82 + ay? + a?) (8), 
and call 7’ the kinetic energy of the particle. Then by in- 
tegration with regard to the time from t to ¢, we get 


[7| == in Xt dt (9). 


But in any short interval of time dt, ¢,d¢ is the increase 
of a,, namely dw, Hence the right side is equal to 
3 yer 
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= [x.ae, taken from the initial to the final position of the 


particle. We call the scalar X,dz, the work done on the 
particle by the force X, in a small displacement da;, Then 


3] Xide, is the total work done by all the forces on the 


particle during the motion. We have therefore the scalar 
relation that in any motion of a particle 


Increase of kinetic energy = Work done on the particle 
(10). 
In a system of particles we may add up the equations (10) 
for all the particles. We now take 7’, the kinetic energy of 
the whole system, as the sum of those of all the particles. 
Then we get 


Pt [eewl-mfisan on 


the double summation implying summation for all pairs of 
particles; the first summation is for the particles producing 
the forces X, and the second for the particles acted on. It 
may happen that the X, are all functions of the coordinates 
alone, and not of their velocities, and that provided the 
initial and final values of the x, are the same the integral is 
the same however the x; may vary in the interval. If so, 
the integral is the difference between the values of a certain 
function U for the initial and final positions of the system; 
we call the system conservative and U the work-function. 


‘Then (11) becomes 
P-L os 


or T — U = constant (13). 


This is the equation of conservation of energy. The quanti- 
ties on both sides are scalars. — U is often denoted by V 
and called the potential energy. 
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If U exists, then in all possible small displacements of 

the system 

TEX,dx, = dU (14), 
where the dx,, 3n in number for n particles, are all in- 
dependent. The coefficient of da, for any one particle is 
=X,, the total force in the direction of x, on that particle. 
Hence the force on a particle in the direction of x, is 
@U/éx,, which is of course a vector. 

We can also take the scalar product of (1) by any set of 
small quantities 5x, whatever that constitute a vector. 
Then 

mE, Sx, = UX ,dx; (15) 
is a scalar relation; but as the 52, are arbitrary we can 
equate their coefficients and regenerate the equation (1). 
If we now integrate with regard to t from , to t, we get 


fimadeydt = [* Ex ,de,at 16 
partons | eer (16). 
The left side, on integration by parts, gives 
A La fed ad 
[made] — f* may $8, dt (17). 


But we can consider the x, + 5x, as coordinates of a par- 
ticle in a motion differing slightly from the actual one; that 
is, at a given value of t, the coordinates are x, + 52, instead 
of 2,. Then 


Sox = 9 (a, + 82) — 4 
= &, + B%, — a, 
= ba, (18), 
since @, + 6%, is simply the varied velocity or rate of change 
of the varied coordinate x, + 5x,, Then 


[im 4, Bat = f mé,32dt 


te 


~ [23 mea) dt +0 (68)? (19). 
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Then for every particle, to the first order, 


[i 0 Gmée) + 2x,b2) ae= [mé,2,]* (20). 


We can now add up these equations for all the particles of 
the system. If U exists we can express the result in the 
form 


fis (0+ U)dt=> [mad] (21), 


or, since the limits are not varied, we can move the 6 out- 
side the integral. If the variation is such that the initial 
and final positions of the system are unaltered, dz, = 0 
when ¢ = 4, or ¢,, and we have 


afars U)dt=0 
to 


to the first order in the variations of the path. This is 
Hamilton’s principle. 


EXAMPLES 


1. If the Cartesian coordinates of every particle of a system are 
known functions of a set of generalized coordinates g,, prove that 


ama,te,~ ff Gi) 5) 
3X by = 3X57 tye 


where the summation convention is also understood on the right 
side, Deduce Lagrange’s equations for the case where the g, are 
all independent. 

2. The equations of motion of a particle are ma, = X,; — kt, 
where k is constant. Prove that 


e a 
QT = — aX, + ap (tnt?) + gy (hee). 


Hrence show that for a system in periodic motion, or in one slowly 
changing its state, on an average over a long time, 27 = — Xz,X;,. 


CHAPTER IV 
DYNAMICS OF RIGID BODIES 


A rigid body is one such that whenever it is displaced the 
distance between any two particles of it is unaltered. Since 
three particles A, B, C are in a straight line if the sum of two 
of the distances AB, BC, CA is equal to the third, it follows 
that straight lines are unaltered by displacements of a rigid 
body. Since when A, B,C are not in a straight line the angles 
of the triangle ABC are determinate functions of the three 
sides, it follows that all angles are unaltered by displace- 
ments of a rigid body. If three lines meeting at a point and 
fixed in the body are mutually perpendicular before dis- 
placement, they are still perpendicular after the displace- 
ment. 

The equations of dynamics in the form mz, = X, are 
true with respect to dynamical (that is, non-rotating and 
unaccelerated) axes. Let any particle Q of a body have 
coordinates 2, with reference to dynamical axes at O. Then 
let the body be displaced in any way, and let the particle 
have the new coordinates x,'. We require the relation 
between x,’ and z,. Consider a particle of the body, P say, 
whose coordinates before and after the displacement are 
a,anda,. Put 

Haat ys =a! + yl (1). 
Then y, and y,’ are the coordinates of Q with respect to 
axes at P parallel to the dynamical axes before and after 
the displacement. Also if we imagine the original axes at 
P to be specified by the particles on them, these particles 
in the new position still specify a set of rectangular axes, 
with respect to which the coordinates of Q are still y,. If 
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the cosine of the angle between the y, axis in its new posi- 
tion and the y; axis in its old one is a4, we have therefore 


UW = GaeYe (2), 

whence ay! = Oy + Oy (% — MH) (3). 
The displacement of Q is 

a! — a + Gate — Ye (4). 


We can prove that there is a straight line of particles in the 
body such that the displacement of Q is the same as that 
of P. For such a particle we should have 
Gn Ye — Yi = (ie — Bix) Ye = 9 (5); 
and these three homogeneous equations in y, are consistent 
provided 
Wa —8a||= (6), 
that is, q—-1 yy dy =0 (7). 
Gq Ay— 1 das 
Oy, Gq yg — I 
Now the determinant || a, || is unity, and each element in 
it is equal to its first minor. If we expand (7) we get 

Wa || = {(@ix@22 — Gazer) + (4e2%s9 — G29 %s2) 

$+ (Aggy — 51%s)} + (Gar + Gee + a)- 1= 0 (8), 
since the terms all cancel. Hence (5) have an infinite 
number of solutions, all proportional, and the points there- 
fore lie on a straight line. 

If we take any two planes through this line, the angle 
between them is the samo after displacement as before, 
and therefore all planes are rotated through the same angle. 
Thus any displacement of a rigid body is equivalent to a 
displacement of a particle of it combined with a rotation 
about a line through that particle. If the angle of rotation 
is 0, we have, by comparing (5) with (52) of Chapter 1, 

bun = Ox — be 
= (hile — 8a) (1 — 008 8) — €pmlm Sin 8 (9), 
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where J, are the direction cosines of the axis of rotation. 
Thus the a, are determined in terms of this axis and the 
angle of rotation. Conversely, if we take the symmetrical 
and antisymmetrical parts, 


(Lh, — 8a) (1 — 008 8) = $ (bax + Das) (10), 
4 (xim — Edm) Im 810 O = 4 (bie — bux) (11), 
But for a given m, with ikm in cyclic order, 
ti toe a? (12). 
Thus Ty sin 0 = — } (bu — byi) (13), 
and 
sin? 0 = § {(b12 — byr)® + (bos — by2)® + (Bsr — bra)*} (14), 


Now suppose that the displacement is small. Then a, 
have nearly their values for zero displacement, that is, 
5,- The direction cosines of the y, and y, axes in their new 
positions are a, and @,,, and thus, if k # l, 


A.A = 0 (15), 
and, if k =1, nA = 1 (16). 


But in (15) for i = k, ay is nearly 1 and ay, small; for i = 1, 
a, is nearly 1 and ay, small; for i not equal to & or 1, both 
@,, and a,, are small. Hence to the first order 


ay, + y= 0 (17). 


Ifk=1, then for i # k or! both terms of (16) are small of the 
second order, and therefore for i= k or 1, ay, =1+a 
second order quantity. Thus to the first order 


Ga = 8a — ba (18), 


where 6,, is an antisymmetrical tensor. 
‘The displacement of Q is 8a, + a4 — y,, where Sa, is 
the displacement of P, and is a first order small quantity. 
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Now let X; be the total force acting on the particle at Q. 
In the displacement the work it does is the scalar 


BW = X, (8a, + day — ys) 
= X, (6a; — bay) 


= X,8a, — Bae (Xiyn + Xeys) — Ww (Xie — Xeys) 
(19). 


When i and & are interchanged, by, and X,y, — X,y; are 
reversed in sign, while X,y, + X,y, is unaltered. Hence 


8W = X80, + 4 (ys Xe— Ye Xy) ba (20). 
Now write YW Xn — WX = Lay (21), 


so that Z,, is an antisymmetrical tensor. The second term 
in (20) is the sum of nine terms, of which three are always 
zero and the others equal in pairs. If we replace by, and 
Ly, by the associated vectors, we have 


BW = X80; + Lin dm (22), 


where the Sa, and 6,, are all independent of one another 
and the same for all particles of the body. Hence if we add 
for all particles of the body 


BW = (EXj) Bay + (ZL) Om (23), 


and the total work done in any given small displacement is 
determinate if we sum up the forces acting by the six 
expressions 2X,, LL,,. Further, the contributions to 
=X,, XL,, from the internal reactions are zero. This follows 
at once if these reactions consist of equal and opposite 
forces between pairs of particles along the line joining 
them, and also has the justification that it leads to correct 
results. Then we may restrict 2X,, UL,, to the contribu- 
tions from the external forces. This is d’Alembert’s Prin- 
ciple. 

In the limiting case of continuous motion, we may con- 
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sider the displacements that take place in a short time 
interval St, and put in the limit 


8a, = u,8t; bg = wy dt (24), 
and call u, the velocity of P and w,, the angular velocity 
of the body. Then the velocity of Q is 

By = Uy — WEY (25). 
Now consider the centroid @, with coordinates % de- 
fined by 


(2m) B, = Uma, (26). 
It is usually assumed without proof that G is fixed in the 
body, though this is not obvious. But suppose that the 
particle at @ with coordinates 2, in the original position of 
the body goes to @’ when the body is displaced, its new 
coordinates are 7,’, given by 


Bl = aj + ay (7 — &) (27), 


and the coordinates of the new centroid @” are 7,”, given by 
(2m) @," = Ima’ (28). 
We have to show that @’ and @” coincide. We have 
(Sm) Z" = Lm fay! + ay (Xe — a%)} (29), 
and therefore 
(Zam) (%4" — F,') = Sm fay’ + ain (@ — a)} 
— (Bm) {ay + ay (7% — a)} 
= Lmay et, — UmayT, 
=A (30). 
Thus the particle originally at @ is displaced to the new 
position of the centroid, and therefore the centroid is fixed 
in the body. 
Now return to the equations of motion of a particle of 


the body 
mii, = X, (31). 
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By addition we form from these the equations 
Lmé, = BX, (32), 

and by cross-multiplication followed by addition for all 
the particles of the body, 

Lm (aay — ty H;) = U (aw X,—%Xi) (38). 
On the right sides of (32) and (33) the contributions from 
the internal reactions are zero, by d’Alembert’s principle. 
Also 


ner: 2 - 
Yi, = OF Ema, = ue (&m)z,= Ma, (84), 


where M = =m (35), 
the total mass of the body. Also if 
Tea = Lan (20a, — yts) (36), 
we can reduce (32) and (33) to 
Mi, = 2X, (37), 
‘ ha = SL (38). 


These are the fundamental equations of rigid dynamics. 
The three independent h,, are expressible in terms of the 
three independent w,, by (36), and we have therefore six 
differential equations for the Z, and w,, which determine 
them in terms of the initial conditions if the external forces 
are given. The motion of the body is therefore determinate. 

The principle of virtual work follows immediately. For 
if a body is initially at rest, 7, and hy, are zero, and the 
condition that they may remain zero is that 2X, and LL. 
shall vanish. But this implies, by (23), that in any small 
displacement of the body the work done by the external 
forces is a small quantity of the second order in the dis- 
placements. Conversely, if there are six independent 
possible small displacements such that the external forces 
do no work, to the first order, in any of them, the co- 
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efficients of the independent Sa, and 6,, of (23) must all 
vanish, But these coefficients are the XX, and ZL,,, and 
therefore the ¢; and hy do not vary with the time. The 
vanishing of the work done by the external forces in six 
independent small displacements is therefore a necessary 
and sufficient condition for equilibrium. 

The equations (38) may be put in another form, If we 
consider any moving origin O', not necessarily fixed in the 
body, with coordinates a,, we can write for the coordinates 
of a particle with respect to O’, 


Y= Be — a (39). 


If any vector associated with the particle, such as its 
velocity, momentum, or acceleration, or the force acting 
on it, has components u,, we may form the antisymmetrical 
tensor y;, — y,u and call it the moment of the vector 
about 0’, From (32) and (33) we can form the equations 


Em (x44, — %¥,) — (a, Lme, — G, Zms,) 
= ¥ (aX, — %X,) — (EX, —a,EX,) (40), 
where only the external forces make any contribution to 
the right side. But by (39) this is equivalent to 
Lm (yiiie — Yui) = B (ys Xe — ye Xe) (41). 
Therefore the sum of the moments of the mass-accelera- 
tion products about any origin is equal to the sum of the 
moments of the external forces about that origin. 
If we denote the moment of momentum, or angular 
momentum, about O’ by hy’, we have 
ha! = Em (Yb, — Yue) (42), 
and 
Bisey a . Sees 
a hg! = Em (yee — Yxts) + Lm (Yite — Yeti) 
(43). 
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‘The second term may be written 
Bim {(v; — 4) by — (ee — Gy) &} = — Lm (dj, — G4) 
=— M (ij, — a3,) (44). 
Thus (41) is equivalent to 
d 
a 
In many important cases the second tensor on the left 
vanishes identically. This is clearly true if O’ is fixed, when 
the time derivatives of the a, are zero; when the centroid 
is fixed; and, if both are moving, if the vectors a, and @, 
are parallel, that is, if the velocity of the moving origin 0’ 
is parallel to that of the centroid. The most important case 
is where the moving origin is identical with the centroid, 
when the last condition holds automatically. These terms 
also disappear if the moving origin is an instantaneous 
centre of rotation always at the same distance from the 
centroid, They vanish for a sphere or circular cylinder 
rolling down an inclined plane, but not for a rolling elliptic 
cylinder. 


Ign! + M (tts, — Gt.) = Z(yXe— 2X1) (45). 


If the moving origin is the centroid, 
ha’ = Em {yi (oe + Yr) — Yu (be + Gd} 
= Lm (Ye — Yeti) (46), 
since Xmy, = Imy, = 0, 


by the definition of the centroid. Thus the angular mo- 
menta about the centroid are expressed completely in 
terms of positions and motions relative to the centroid, and 
the formulae for them have the same form as those for the 
total angular momenta with reference to a fixed origin. 
By (25) Ke = — mY (47), 
m being here a dummy suffix. Substituting in (46) we have 
Tia! = — EM (em YY — ime Ym) 
= Bemeoim — BinOxm (48), 
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where B,,, is the symmetrical tensor defined by 
Bom = EMY Yn (49), 
and depends only on the masses and positions of the par- 
ticles of the body. In the three-dimensional case, with i 
and & unequal, m must be equal in turn to i, k, and the 
other value. Then, for instance, 
Iya! = (Buon + Brey + By) 
— (Buon + Byytog, + Byycros) 
= (Ba + By) Oy — Bysry — Byyeogy (50). 
If we replace the antisymmetrical tensors by the associated 
vectors this takes the form 
Ig! = (By + By) os — Byer, — Bryeog 
= (By + Buy + Beg) oy — (By, + Byywg + Byyers) (51), 
and in general 
Ion! = Agen, (52), 
where A,,, is the symmetrical tensor given by 
Ain = Budim — Bim 
= my," 8m — SY Ym (58). 
Tt evidently corresponds to the system of moments and 
products of inertia given in ordinary treatises on dynamics; 
Ay, Ay, Ag are the ordinary moments of inertia A, B, C, 
but 4,3, Ay, Aj, are equal and opposite to the ordinary 
products of inertia F, G, H. 

It should be noticed that this reduction is characteristic 
of three dimensions; in a higher number of dimensions 
there is no analogous simplification of the form (48). 

The equations of motion of a rigid body then take the 
form 

Mi, = 3X, (54), 
d 


Gpha’ =ELa'; or Shy’ =EL_’ (68), 
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where L,,’ is the moment of the external forces about the 
centroid. 

It may happen that one point of the body is fixed. In 
that case we may take the origin at that point. Then in 
addition to the known external forces there is a reaction 
at the origin which can if required be found from (37); but 
the reaction has no moment about the origin and the 
motion is given by (38). But in this case 

£5 = — WimFZm (56), 


and we find, by a process analogous to the last one, 


Tie = Dim im — Dim em (57), 

where Din = TMX j%—q (58). 
In terms of the associated vector, 

Tin = Comeng (59), 

where Cim = (TMX2) F4m — IMT; Ty, (60). 


The C,,, correspond 0 the moments and products of inertia 
about the origin. They can be expressed in terms of those 
about the centroid; for 


Com = {Em (Ty + Yu)"} Bien — Em (e+ Ys) (Zn + Yon) 
= ([mZ,") 84m — UME, 7, 
+ Am (61), 


the terms linear in the y’s vanishing by the definition of the 
centroid. 
‘The relevant equations of motion then take the form 


d 
qin = Dn (62), 
where the Z,, are the moments of the external forces about 


the origin, the reaction at the origin making no contri- 
bution. 
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The kinetic energy of the body is given by 
27 = Ime? 

= Lm (& — Hmm)? 

= Lie + TMP YmOnom 

= Mi? + Bymoneim (63). 
In three dimensions w .«;,, is zero unless ¢ is different from 
both & and m. If i= 1, k = 2, m= 3, or if i= 1, k= 3, 
mM = 2, wets = — WH. If i=l, k=m=2, wy? = w,%, 
and has coefficient B,.. Thus w,* enters with a coefficient 


By + By, or Ag. But w,w; has a coefficient — 2B,, or 
Ag + Aq. Thus in all 


27 = Mi* + Ayo (64). 
When a point of the body is kept fixed, 
27 = Opn Opn, (65). 


We notice that the linear and angular momenta take the 
forms 


4 Peabo 
Mi, = a h{= Jo, constant), 
or 
hy= ose (66) 


for origin at a point of the body held fixed. 

In one respect the forms (55) and (62) are inconvenient. 
They involve the tensors A,,, and C,,,, which depend on 
the x, and y, and therefore in general change as the body 
rotates. It is more convenient to use such axes that the 
relevant tensor in the particular problem is constant. To 
achieve this the axes must rotate, and then are no longer 
dynamical axes. Suppose then that we have a set of axes 
a,’ rotating in any way, and that their direction cosines 
with respect to the dynamical axes x; are a,;. All the usual 
tensor relations hold for transformations from the 2,’ 
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system to the x, system and vice versa; but the a,;, while 
specified for each instant, are now functions of the time. 
If u, is any vector, we have 


Uy = Ay Uys Uy = Oss (67). 


Then a is another vector; and 


= ats + asa (68). 
The component of this vector in the «,' direction is 
Ay Qyy ay + yu; ee (69). 
The first term of this is 
bye = Set (70) 


Also ee is the x, velocity of a point at unit distance along 


the «,' axis, that is, a point with coordinates a,, (i = 1, 2, 3) 
with reference to the dynamical axes. This velocity is 
— 04,3, where 0,, is the antisymmetrical tensor expressing 
the rotation of the rigid frame consisting of the moving 
axes. Then 


da, 
1 Mays 
Carty ag ei js Oe Aey 


= — Oy us (71), 
and the required component is 
ses — 6y'u, (72). 


If we use instead of 0,’ the associated vector, the three 
components become 


(ig! = t40,! + ty’ By, tig? — 14" Oy! + 4", 
ti — 14 0,' + u4'0,!) (78). 
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If for instance the components of a displacement along the 
2,’ axes are u,’, the formula (73) gives the components of 
the velocity along these axes; if u,’ are the components of 
velocity, (73) gives the components of acceleration; if u,' 
are the components of angular momentum, (73) gives the 
components of their rates of change with reference to 
dynamical axes, and these are equal to the components of 
the moments of the forces acting. 


EXAMPLE 


Prove that 
hig! = — SMI Cimy mks Ys 


= Aina. 


CHAPTER V 


EQUIVALENCE OF SYSTEMS 
OF FORCES 


An external force X, acting at a point 2, of a rigid body 
produces dynamical effects summed up in the vector Xs 
and the antisymmetrical tensor z,X,—2,X,. A force has 
a line of action; that is, if we take its resultant R given by 


R=X? (o)) 
we can define a direction 1, by 
X, = Rl, (2). 


By convention R is always taken positive. If a foree X, 
acts at the point x, + rl,, where r is arbitrary, we have 
(ay + ly) Xp — (ate + rh) Xp = UX, — YX, (3). 

Thus the dynamical effects are the same if the force X, 
acts at any point of a line through 2, with direction cosines 
proportional to X,. This line is called the line of action of 
the force, and the force can be said to act along it. 

If a force has magnitude R and acts at z, in the direction 
1,, we have 


X,= Rl, (4), 
Ly = RB (ihe — %li) (5), 
or Ln = Bln’ (6), 


in the notation of the coordinates of a line. Thus the X, 
and L,, are the products of the resultant into the six 
coordinates of the line of action. 

The moment of the force X, about a point a, is 


(@_— @;) Xp — (@e — %) Xp = RB (lx’ — ale + ali) (7). 
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If we take the associated vector, its component in the 
direction m,, is 
R (Malm! — rm Mm Aly) = B (ml! + lm,') (8), 


where m,’ are the other coordinates of the line through a, 
in the direction m,. We notice that 

ml! + lm = dsin 0 (9), 
where @ is the angle between the two lines and d the length 
of their common perpendicular. We may call Rd sin @ the 
moment of the force about the line (m,, m,’). 

By d’Alembert’s principle, the motion of a body is un- 
altered if to the forces acting on it we add two equal and 
opposite forces acting at the same point, or, by (3), along 
the same line of action. 

Now consider a pair of equal and opposite forces X, and 
— X;, acting at points a; and b,;. They clearly make no 
contribution to £X,. Their contribution to LL, is 


(a, Xp — aX) — (0, Xe — bX) 
= (a — b) Xp — (@ — by) Xe (10). 


Since a, and 6, are equally affected by any motion of the 
origin, the contribution of such a pair of forces applied to 
definite particles to both 2X, and EL, is independent of 
the position of the origin. Such a pair is called a couple, 
and its contribution to ZL, is called the moment of the 
couple. 

If the vectors a, — b, and X, are both perpendicular to 
a line with direction cosines n;, the components of EL,, are 
the components of a vector along this line; while the forces 
act in the same plane perpendicular to the line. This plane 
is called the plane of the couple, and the line an axis of the 
couple. Evidently equal couples in parallel planes are 
equivalent. The magnitude of this vector is Rd, where 
R? = X/ and d is the perpendicular distance of a, from a 
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line through 6, parallel to X,; and its components can 
accordingly be written Rdn,,. 

Any system of forces is equivalent to a force at an 
arbitrary point a, together with a couple. For with each 
X, acting at x, we can associate a pair of forces + X, at a;. 
Then the system is equivalent to SX, at a, together with 
a set of couples whose total moment is 


My = E {(eq — 4) Xe — (ey — A) XH (11). 


But this is equivalent to a single couple; for we have only 
to make n,, proportional to the M,, and Rd equal to their 
resultant. 

If Ly are the moments about the origin, 


Diy = Be, Xy — %Xy) (12), 

and if for brevity we replace 2X, by simply X,, 
Ma = Liye — (a X_ — GX) (13), 
or My = Lin — tem Xp (14). 


Evidently X,* is a scalar and independent of a,. Also 
Xn Mn = XnLin — em Xe Xm (15). 
The first term is a scalar and independent of a,. The second 
is identically zero. Hence X,* and X,M, are scalar in- 
variants. 
The system is equivalent to a force X, and a couple M, 
at a,;. These vectors are parallel if 


Lim = €tem% Xz = PXm (16), 


where p is a scalar length called the pitch. These give three 
linear relations between the three a, and the pitch, and we 
therefore expect a single infinity of solutions. But if we 
take the scalar product of (16) by X,, we have 


Xm Lm = PXm* (7), 
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so that p is the same for all admissible values of a,. With 
this value of p, (16) represents three planes with a line in 
common. If we change a, to a, + ¢X,, where a is a scalar, 
the left side of (16) is increased by 

=n OX Xe = 10 (18), 
and therefore if a, is one point on the line, all points on the 
line through a, parallel to X, satisfy the conditions. This 
line is the central axis of the system; the system is equiva- 
lent to a force R along the central axis and a couple @ 
about it. This expresses the system as a wrench. Tf we take 
the central axis as one of the coordinate axes, we have, 
since X,* and X,M, are scalars, 

R=Xe (19), 
GR = X,M, = pX? (20), 
and therefore G, R, and p are determined. 
The system can also be reduced to a couple parallel to 
a preassigned plane together with a force. For if 8 is the 
couple, and n, are the direction cosines of the normal to the 
plane, and if the force acts through a,, we have, for the 
moments about a,, 


Mu = Dei €omiXn (21), 
and also My = Sin (22). 


We have three equations to determine the a, and 8. Again 
there are a single infinity of solutions. But if we take the 
scalar product by X,,, we get 
Sty Xm = Lin Xm = im Xp Xq 
=GR (23), 
so that § is determined provided z,, X,, is not zero, that is, 


provided the resultant force is not parallel to the plane. 
‘Then the equations 


Lin — €ixm% Xp = Stim (24) 
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determine a line parallel to X;, which is the line of action 
of the force. 

If 1, are the direction cosines of a line through a;, the 
moment of the system about this line is 


Moe ear Xk, (25) 
= In Lim + hy’ Xe (26), 
where J,’ are the other coordinates of the line. If this 


moment vanishes the line is called a null line of the system. 
If 6, is another point on it, 

Ly (Om = mn) = € cx Xx (Om — Om) (27), 
which shows that , lies in a definite plane through a,. All 
null lines through a point therefore lie in one plane. This 
plane is called the null plane of the point. 

All null lines in a plane pass through a point. For the 
system can in general be reduced to a couple in the plane 
and a force whose line of action intersects the plane in one 
point. Then the system has no moment about any line in 
the plane through this point, which is the null point of the 
plane. 

Any system is equivalent to two forces, one of which can 
be made to act along a given line. For let the lines of 
action pass through a, and 6,, and have direction cosines 1, 
and m,, and let the magnitudes of the forces be S and 7’. 
‘Then we have six equations, 


X, = Sl, + Pm, (28), 
Lin = €ixm (4/Sh, + b,Tm,) 
= Sly! + Trt! (29). 


The coordinates of the first line being given, these are six 
linear equations to determine the six coordinates of the 
second line and § and 7’. But we have also 


m2 =1; mm =0 (30), 
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so that we have in general just enough equations. Two such 
lines are conjugate lines. Clearly any line intersecting two 


conjugate lines is a null line. 
It can be shown easily that 
S(X,li + Ljl;) = LX; (31), 
T? = X?— 281,X,+ S* (32), 


whence (28) and (29) determine the coordinates of the 
second line explicitly. 


EXAMPLES 


1. A system of forces is reduced to a force at P together with 
a couple; P is chosen so that the couple is parallel to a given 
plane. Show that the locus of P is a straight line parallel to the 
central axis. 


2. Show that 
a, = SiemXily 
im 
“ 
is a point on the central axis. 
8. Two systems of forces are given by (X,, L,),(¥,;, K,). Show 
that X,K,+ Y,L, is invariant. 


CHAPTER VI 


CONTINUOUS SYSTEMS 


In problems involving volume and surface integrals we find 
it convenient to denote elements of volume and surface by 
dr and d8 respectively; both are always taken positive. 
They are of course scalars. The direction cosines of the 
normal to an element of surface are usually denoted by J,; 
in most cases the normal is drawn outwards from the region 
under consideration. 


Green’s Lemma takes the form 
(e-ffinae oy 
and we have the corollary, if 
u, = OV /dx, Q), 
where V is a scalar, 
Stereo 
where Va ae mat mat m= (4), 
and @/én denotes differentiation along the outward normal. 
Stokes’s Theorem takes the form 
Jmiee= [ff -S2) + 4G 3 = 
+h (g- mf as 
-| | Lem Otlyg/2ty AS (5). 


The integral on the left is round a closed contour C. On 
the right J, is the direction cosine of the normal to any 
element dS of a surface S that fills up the contour; the 
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integral is over S; and the sense of the normal is such that 
if the contour is described in the positive sense about any 
axis z,, 1, is taken positive when the normal is in the direc- 
tion of x, increasing. 

The gravitational potential of a distribution of matter is 
given by 


v-sa2 6), 


where m is an element of mass, r is its distance from the 
point where V is to be found, and f is a constant equal to 
6-66 x 10-8 when m and r are measured in grams and centi- 
metresand the time in seconds. When themassis distributed 
continuously over a surface or through a volume, m must be 
replaced by od or pdr respectively, where o and p are called 
the surface and volume densities. The work in displacing 
a mass m’ through a small distance is m/dV. 

The electrostatic potential of a system of point charges 
is given by 

V=sur (1), 


where e is a typical element of charge. If ¢ is in electro- 
static units of charge and rin centimetres, fis + 1. The work 
ina small displacement of a charge e’ in the field is— e’dV. 

The usual relations follow, that in free space in both 
cases 


vey =0 (8), 
and in space occupied by matter of finite density 
V°V = — 4nfp (9). 
Also we have Gauss’s Theorem 


Kz oF d§ = — 47fX'm or — 47fd’e (10), 


where the summation is over all the masses or charges 
within the closed surface S. In crossing a surface where 
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there is a finite surface density 0, V/én has a finite dis- 
continuity — 4zfc. 
All these relations are scalar in form. For a gravitational 
field the force on a small particle of mass m’ is 
m'X, = m' dV /dx, (11), 
while for an electrostatic field the force on a small charge 
eis 
eX, = — e' dV /da, (12). 
The vectors X, are called the intensities of gravitational 
and electric force respectively. The analogy in form be- 
tween (6) and (7) is constructed for mathematical con- 
venience; the difference in sign in (11) and (12) embodies 
the physical difference that whereas two positive masses 
attract, two positive charges repel. 
In a gravitating system we may construct a work- 


function 
MyM 


W=,2f (13), 


Toa 

where the summation is over all pairs of particles m,, m4, 
and r,, is the distance between m, and m,. The case where 
p= is excluded. Then the force on the pth particle is 
@W /(ea,)p- The potential at m, due to the other particles is 
V, = 2'fmq/rnq, the accent indicating that the case where 
q= p is excluded from the summation; and the function 
Im, V, = 2W, since each pair of particles is counted twice 
in this double summation. The function 32m, V,, therefore 
plays the part of a work-function. Similarly in electro- 
statics the function }Ze,V, plays the part of the work- 
function with its sign reversed, thatis, of a potential energy. 
Those results may be generalized to the case of continuous 
distributions; thus we can replace these functions by 


w= aufffevar Kd iffovas (14), 
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where the first integral is through all regions where p is 
finite and the second over all surfaces where a is finite; the 
first integral may therefore be taken through all space, 
excluding the surfaces where there is a surface density. It 
can be shown that with proper precautions about the de- 
finition of V the restriction p # q gives no trouble provided 
V is everywhere finite, 

Now if we consider the integral through all space except 
thin laminae surrounding the surfaces where g is finite, 


afffevar =- aap JIN? aaa 
~~ alll &(7)-G} 


1 ov 1 avy? 
~ sap\] 5% * sap Ie) @ 09) 
where the dn in the first integral is out from the region of 
integration and therefore towards the surface where o is 
finite. On the two sides of such a surface the values of V 
differ by an indefinitely small amount, and for the two 
av 


sides together an = 4nfodS, by Gauss’s Theorem. 


Hence the first integral in (15) is equal to — aI oVds 
taken over the surface and cancels the second term in (14). 


pe aylff(e 
= set I) Ge, (16) 


through all space. In consequence of this form we may say 
that the gravitational work-function, or the electrostatic 
potential energy, is R*/8zf per unit volume, where F is the 
resultant of the appropriate intensity vector. 

When the properties of the medium vary from place to 
place, V is no longer of the form Sfe/r, and V?V is no longer 
zero. But a potential still exists; if a small charge e’ is 
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moved from a point P to a point Q, the work done is still 
the same whatever the route, and may be denoted by 
é (Vp— Vo). The treatment is suggested by the fact that 
two similar condensers with the plates at the same po- 
tentials, but with air between the plates of one and another 
material between those of the other, have charges in a ratio 
K depending only on the media. @V/én on the outside of 
the condenser being small compared with its value be- 
tween the plates, we infer that the charge per unit area for 
the same distribution of V is related to the discontinuity in 
KoV/én, where K depends on the material. This suggests 
in turn that Gauss’s Theorem must be replaced by 


\ja Zas- — dafBe (17), 


where the summation is for all charges inside S. Then 
applying this to the two sides of any surface we have 


[« rl espe (18), 


and applying it to a region with a finite volume density we 


have 
-saifffote= ff (Ze 0 


and as this must hold for all such regions, 


é (K ar) so Bere (20). 


62, 


These equations are all homogeneous in V, p, a; hence 
the potential due to any set of charges is proportional to 
the charges if all are altered in the same ratio. Using this 
principle we can show by the usual method that the energy 
of a distribution is 


y2Vve=4 [fear + 4 [fovas (al). 
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The first integral, applied to all space except thin regions to 
cut off surfaces where there are surface densities, gives 


~ eal) "i (Ka) 


“i J fee, ( 


ve) -& (er) ‘} a 
~~ aap JEP 35 29+ apf] Ge) & 
The surface integral cancels the integral sfforas 3 and 


hence 
Ws lllt Ge) 


so that the energy can be considered equal to KR?/8xf per 
unit volume. 

Magnetism may be treated similarly, starting from the 
assumption of volume and surface distributions of mag- 
netic pole strength, subject to the condition that the total 
pole strength in any solid is zero; or we may regard the 
ultimate magnetic unit as the doublet, which explains the 
need for the restriction involved in the former method of 
treatment, The potential at x, due to a doublet of strength 
M at the origin with its axis in the direction A, is 

V = yMA,x,/r8 (24), 
and if the doublet strength per unit volume in a solid is I 
in the direction , we can introduce the intensity of mag- 
netization at ¢;, the vector A,= JA,, and say that the 


potential at a, is 
V= fall [4s f+ ap (25), 


where A, corresponds to the point £, and dr = dé, dé,dé,; 
y is a constant. The magnetic force in free space is 


eee he (26). 
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We may write 


var {ffacze (7) * 
ptt [Mee (4) 2564 
nrfitan-rffite e 


The potential is therefore equivalent to that due to a dis- 
tribution of magnetic poles 1,4; per unit area over the 
boundary and — 04,/0é, per unit volume through the in- 
terior. 

Within a solid special treatment is needed. To define V 
or the force at z,, when 2, is within a solid, we must imagine 
a small cavity made about «,, the intensity of magnetiza- 
tion everywhere remaining as before, and consider V and 
«, within it; then the values of V and a, at «, are defined to 
be the limits of those in the cavity when the dimensions of 
the cavity become indefinitely small. This process leads to 
little difficulty in gravitational and electrostatic problems, 
but in magnetism the limit of the force is found to depend 
on the shape and orientation of the cavity. The force in 
the cavity can be written 


<j== 5 (28), 


where V is given by (27); in the first integral the normal is 
inwards towards the cavity. The contributions to X, from 
the volume integral and the outer boundary are of the same 
form as for gravitation, and give no trouble. If the cavity 
is a cylinder in the direction of the intensity of magnetiza- 
tion, 1,4; is zero over the sides and equal to J, the resultant 
intensity of magnetization, on the ends. Such a surface 
density over the ends in the limit contributes nothing to V; 
if the radius of the cylinder is small compared to its length 
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it also contributes nothing to X,; but if the cylinder is 
of disc-like form it contributes 4zyA,J to X,;. Hence if we 
take V for the complete body and define a, by 
av 
os a (29), 

a, is the value taken by X, in a thin cylindrical cavity 
parallel to the intensity of magnetization. The force in a 
flat cylindrical cavity with its generators in this direction is 


a, = a + 4ryAI = a + 4ryA, (30). 


Evidently «, and a; are both vectors; the former is called 
the magnetic force and the latter the magnetic induction. 
The theory of susceptibility and permeability may then be 
developed as usual. Also (27) shows that V is continuous 
across a boundary; but J,a;= 2V/dn has a discontinuity 
— 4myl,A,;; whence 1,a, is continuous across a boundary. 

The mutual potential energy of two doublets M and M’ 
at x, and z,’, oriented in directions A, and A,’, is 


ev 


W= WM 5 (31), 


where V is the potential at x,’ due to the magnet at «,; this 
gives 


W = ,MMa! e (a 


= yar {Ma 


ae (cos ¢ — 3 cos 8 cos 6’) (33), 


where ¢ is the angle between the axes of the magnets and 
6 and 6’ are the angles made by the axes with the line 
joining the centres. 
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If the second magnet is turned through a small angle 5p 
about a line with direction cosines ”,, 


BAY = eam Mn Am’ Sih (34), 
by (49) of Chapter m1. Hence, by (32), 


we, 
aw = ay 
= patar’ 8 La, — 3.008 0% = * cant Ae! 


so that the couple about a line parallel to the axis of x, is 


MM’ rae 
My = cam dn {a- 3 cos oma (35). 


Hydrostatics and Classical Hydrodynamics. The internal 
reaction in a fluid across an element of surface dS is a 
pressure pd normal to that surface. If the density is p, 
the bodily force per unit mass X,, and the velocity of the 
fluid at x, is u,, the acceleration of the fluid is found, by 
considering a small parallelepiped, to be given by 


du t) 
Pat = — aut eX (). 


If u, is given in the Eulerian way as a function of the 
coordinates x, and the time ¢, the operator d/dt, giving the 
rate of change of any element associated with a given 
particle of the fluid, is equivalent to 

d_2@ C) 
an-at om, (2). 


If we consider the circulation Q around any closed 
circuit C in the fluid and moving with it, defined by 


a (3), 
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a ft + fe y(t) 
-- [12 aa,4 | Kido + [ udu (0. 
e co 


The last integral is [}7,2], which always vanishes because 
when we move round the contour we come back to the 
same point, where the velocity has its original value. Also 
if X, is the gradient of a single-valued potential, as when 
the bodily forces are due to gravity (the commonest case), 
J Xda, is the change of this potential round the contour; and 
is zero. Again, if p is a function of p only, as in an incom- 
pressible liquid or a gas at uniform temperature, the first 
integral vanishes and 


aQ 
7° (5). 


If then Q is ever zero around a circuit it remains so per- 
manently. This is true if the fluid is initially at rest and is 
set in motion by solids moving in it, and in various other 
cases of importance. But the vanishing of Q for all circuits 
is the condition for the existence of a velocity potential 4 
such that 


a. 
uy = ie (6). 
In this case we can rewrite the equations of motion in 
the form 


Oe ge — ae + Xe , 
and multiplying by dz, and adding we have 


ta) Ou, 1 
2 ae, + uagitdy=— 3 pede Kido (8), 
‘ dup 384 _ Oy 
rele Qa, OajOx, Oa (). 
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This shows that for all contemporaneous variations 


ti dj 
ee + dues — [ee U + constant (10), 
where Faye (11); 
i 


u,2 is the square of the resultant velocity q. The constant of 
integration is not necessarily the same at all instants and 
therefore may be a function of the time. Hence we have 
the Bernoulli integral 


4+ i--[P+0+Fe (12). 


The rate of change of mass within a given small paral- 
lelepiped dz,dx,dx, is equal and opposite to the rate of 


outflow; hence 4, a 
BOT = ay, (Od) dr (18), 
and we have the equation of continuity 
a a 
ieee AG (14), 
dp Ou, 
or ime Bis (15). 


Vectors with given Divergence and Curl. We sometimes 
have to find a vector wu, such that 


Ou, 

a, A Q), 
7 Ou, 
il on bamihied (2); 


where A is a given scalar and w,, a given antisymmetrical 
tensor. We want particular integrals of these equations, 


Evidently if Ue = a (3); 
a 
where ¢ is any scalar, 
V4p=A (4), 
Ou, — Ouse _ 
te, come (5), 
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and a solution of (4) is 
1 A 
$--z] [Fe (0), 
where ¢ is to be evaluated at x; and A at €;; 
dr = dé, dé, déy (7), 


and 7 is the distance from ¢; to a. 
_ OF, OF, _ OF y 
If uy Ox, = Bin = Cam Oty (8), 


where F, is a vector such that 


ar, _ 
one ), 
Ou, OF, 
we have Be, = “8 Bebe, 0 (10). 
Then also 
wx Oty _ Oly On 
Om oa, Oty“ Oa, 
C) ar, 
= Stem Be, Expy bx, 
ms oe 
Stim foes Oo Oar, 
oF, 
= (BuBne— Biedm) sea 
Berth iypaamng 
00, O0im "02, 0%y 
0, OB lm 
Om, 0%; Ou,8 
=-V'Fa (11) 
by (9). 
' 1 ([[em 
Thus if Py = qe {f[S2ar (12), 
we shall have 


Seg Ms Be eh ao (13), 
t 


60 CONTINUOUS SYSTEMS 
provided (9) is satisfied. But 


te Ellmd.C)* 


1 Cal 
--glllmaG)* 0. 
since w, is a function of é, alone and r a function of 
x,— €,. Applying Green’s Theorem to all space except a 
small sphere about x,, we get 
2Fm 1 yim ({ mag + 2 tim [[{ 12m 
aoe ~~ dq Lim [| Pas + gta [7 Gee 
(15), 


since F,, has the form of a gravitation potential and 
aF,,/ex, that of a gravitational force. But the first integral 
vanishes in the limit when the sphere becomes very small, 
and the integrand in the second is zero provided the com- 
ponents u, exist, for 

Bum 2 (Atte Ou) _ Ome _ 

Om fa (aa ~ Bey) ~ SO (Bet, 
If then we are given a scalar A and a vector w,, such that 
its divergence is zero, a solution of (1) and (2) is 

a6 , OF, Fy 
“iT det By Be (17), 

where ¢ and F,,, are given by (6) and (12) and ikm are 
in cyclic order. 

This analysis has two practical applications. In hydro- 
dynamics wu, is the velocity and w,, is twice the vorticity, 
denoted by 2£,, in Chapter 1x. Here the divergence of the 
velocity and the vorticity may be given through all space, 
and the velocity (17) satisfies the conditions. If the actual 
velocity is v;, we may put 


(16). 


=U + (18), 
Buu’ _, Sux! _ Btu’ 
and then BOs Be ee =O (19). 
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Thus u,' is the gradient of a scalar ¢’ satisfying Laplace’s 
equation. There is no such scalar that makes the velocity 
finite everywhere, including at an infinite distance, except 
a linear function of the coordinates. Hence u,’ is the same 
everywhere. 

If there are solid boundaries or free surfaces at a finite 
distance u,’ may not be constant. If the region where 
vorticity is present does not extend to a boundary, there is 
no contribution to (14) from points outside this region; we 
therefore take (14) through a region large enough to contain 
the whole of the vorticity. Then the surface integral in (15) 
must also be taken over the boundary of this region, but 
still vanishes, and (17) is still a solution, But (17) may not 
satisfy the boundary conditions, and then we must add to 
u, an irrotational solution chosen to make the whole velocity 
satisfy them. 

In electromagnetism u, may be the magnetic intensity 
and w,, the electric current across unit surface in a plane 
of x, constant. 

In many cases u, has no curl outside a limited region of 
very small cross-section. This is often true in the motion of 
a real fluid, when the region may be called a vortex fila~ 
ment, and in magnetism, when the region is a wire carrying 
an electric current. The former statement may be expressed 
also by saying that the motion is irrotational outside the 
vortex filament; the latter says that magnetic forces due to 
electric currents have a potential. In either case the integral 


on J ugderg (20), 
o 


taken around a closed circuit, is zero if the circuit can be 
filled up by a surface not cutting the filament or the wire, 
and has a constant value for all circuits that cannot be so 
filled up. The two conditions are mutually exclusive, and 
therefore the critical region itself must be a closed circuit. 
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In each case A vanishes. Then ¢ = 0, and the component 


F,, is given by ¥ 
F,= El/ee (21) 


taken through the critical region. But if we consider the 
contribution from an element between two planes sepa- 
rated by dé,,, and call the element of surface in a plane 
parallel to these dS, we have 


dr = dé,dS (22), 
| fond = [uedé, = 2 (23), 
where the line integral is taken around the boundary of 
the filament. Hence A 


ra=Z 4), 


the integral being taken around the length of the filament. 
Also 


ig 2/9 


(25), 


where ds is an element of length of the filament and J, a 
direction cosine. 

In hydrodynamics Q is the circulation around the fila- 
ment. In electromagnetism the unit current is such that 
if it flows in a circle of radius 1 cm. it produces magnetic 
force 27 at the centre. If we take the circle to be in the plane 
of 2, constant, with its centre at the origin, we have 


ds = d0; r=1; 4% =%=2%,=0 
&=c0s 0, & sin, & = 0; h=—sind, k=cos 6, =0 
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and 
Jam £42 ds = 0 for i = lor2 
and = 2 for i = 3. 
Also for unit current u;, the magnetic force, = 27. 
Hence, in this case, by (25), 
Q= 4z, 
and in general, if the current is , Q = 47J, and 


spol eon ee ds (26) 


yf J rag: 7 den (27). 


This may be transformed by Stokes’s Theorem into an 
integral over a surface with the wire as its boundary; thus 


t= 1 | tncne ee én 87 a8 


~~ 1 [ftmcanean sag, (—) #8 
=~ 1 [tm @mBse ~ 89m) 52-5 (3) @8 
--1 se . cane 


[| 
ay [Im 2 5 z¢ (28). 


But ~ ba 5¢ae (2) is is the magnetic force at x, due to 


a doublet of unit strength at ¢, with its axis in the direction 
1, Hence the force at any point is equivalent to that due 
to a distribution of doublets over the closing surface, with 
intensity I per unit area and directed normally to the sur- 
face. Such a distribution constitutes a magnetic shell of 
strength J. 
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It is customary to assert as a fundamental postulate that 
the magnetic force due to an electric current is equivalent 
to that of a magnetic shell, but I think this course un- 
desirable. The magnetic shell does not exist in nature and 
direct experimental test is therefore impossible. Further, 
though it gives the same force, it does not give the same 
potential; the magnetic potential due to a magnetic shell 
is a single-valued function with a finite discontinuity at the 
shell, while that due to an electric current is a cyclic 
function with no discontinuity except at the wire itself. 

Mutual Bnergy of Electric Circuit and Magnetic Field. 
A magnetic pole m at x, is under a force mu; but 


muy = ~ mt 2. [| ba ge (7) 48 (29) 


and — tn 3e-(?) is the magnetic force normal to the 
‘m 

surface at ¢,, due to the pole. Hence — i = Pa 6) dS 

is the total flux across the closing surface, or through the 

circuit, of the magnetic force due to the pole. Also if V is 

the magnetic potential due to the current, 


“oe (30), 
and therefore the correct forces are given by taking 
y=-1{[ was (31), 


where N is the normal magnetic force due to a unit mag- 
netic pole. This can now be generalized; we say that in 
general the mutual potential W of a current and a set of 
magnetic poles is 

‘ast i Nas (32), 
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where N is now the normal magnetic force due to the whole 

of the rest of the field. This form may be further extended 

to express the mutual energy of two electric currents. 
From this equation we can infer, as is done in the standard 


works of Jeans and Livens, that if I | dS varies with the 


time, the variation generates an induced £.M.F. — o If Nds 
in the circuit. y 

If we have two circuits carrying currents J and J, their 
mutual influence is expressed by the statement that their 
mutual energy is of the form (32), where N is taken to be 
the force near £, due to the current J; points on the latter 
circuit can be taken to be given by x,. Then 

W=1JM (33), 


where M= J Lud (34), 


and u, is the magnetic force at , due to unit current in the 
second circuit. But if F, is the vector potential due to such 
a current, 


bo Ste an $8 (35), 
oF, 
M= i canls 9g" aS (36) 
= | Faden (87), 
taken around the circuit. But 

F,~ [% (98), 

taken around the second circuit; and therefore 
a 


taken around both circuits. This gives the required form 
for the coefficient of mutual induction of two circuits. 


CHAPTER VII 
ISOTROPIC TENSORS 


A tensor is called isotropic if its components retain the 
same values however the axes are rotated. We have 
already obtained three examples, namely 5,,, €im, and 
ins €mpe- 

‘There are no isotropic tensors of the first order. For if 
1, was such a tensor, let us give the axes a small rotation 
expressed by the antisymmetrical tensor c,. Then in the 
new system 

tg! = (Bi — C4) Us = Us — Cas My @), 
and this can be equal to x, only if 
CyyUy = 0 (2) 
for all admissible values of the ¢,;. Thus 
Cut + Crate + Crate = 0 
Cy Uy + Cy9 ly + Orgy = 0 (3). 
Carthy + Coote + Canty = 0 
But cy = Cy = C33 = 0, while c,2, Cys, Cy are independent and 
equal and opposite to the components obtained by inter- 
changing suffixes. Hence (3) can be satisfied only if 
thy = tly = ty = 0 (4), 
and therefore there is no isotropic tensor of the first order 
other than zero. 
If uy, is an isotropic tensor of the second order, 
W' ge = (Big — Cis) (Bur — Cnr) yr 
= Ue = Cy Burts — Cer BisUpn 
= Uy — Celine — Orrin () 
to the first order, for all values of i and k. Hence 
Castle + Ces ttsy = O (2). 
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If i and & are unequal, take i = 1, k = 2. Since ¢,.=cy»=0 
we have 
Cratlen + Crs ae + Carty + Costs = 0 (3), 
and therefore 
Ugg = thy = OF Uy = Uae (4). 
By symmetry w,, is therefore 0 if i 4 k, while a, = Uy, = Ugg. 
If i and k are both 1, we have 
Cratln + Crate + Crete + Crstas = 0 (5), 
which is satisfied since every term vanishes. Hence the 
only isotropic tensor of order 2 is a scalar multiple of 5,,. 
If Ujxm is an isotropic tensor of the third order, 
U' tem = (8e5 — C43) (8x2 — Cx2) (Bun — Cran) Msn (1) 
and therefore, for all values of i, k, m, 
Cs Usp + Ces Vis + Cos Ving = O (2), 
Take i= k= 1. Then 


Crore + Cis Usim + C2%em + Cis tism 


+ Catan + Oma trie + Cmts = 0 (3). 
Now put m = 2 so that cy. = 0. Then 


Une + Uae = Man 
Us + tha = Ue (4). 
Uys = 0 
From the last equation, and by symmetry, tj_ = 0 if two 
of i, k, m are equal and the third unequal. Then by the 
first, Um is also zero if all of i, k, m are equal; and the 
second shows that 
Uae = — Unim> 
If in (3) we put m = 1, every term vanishes, so that (3) 
holds. 
Now in (2) if i, &, m are all different, uw», is zero unless 
j =i, and then c,, = 0. Hence (2) holds. It follows that the 
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only isotropic tensors of order 3 are scalar multiples 
of €ixm+ 

If Uamp is an isotropic tensor of order 4, we have, 
similarly, 

CesUsnmng + CnsUcimp + Cmsttarse + Cnstiems = O (1). 
There are only three possible values for i, k, m, p, and 
therefore at least two of them must be equal. We may con- 
sider separately the cases where (a) two are equal and the 
other two unequal, (b) three equal, (c) two equal and the 
other two equal, (d) all four equal. 

In case (a), take i= k = 1, m= 2, p= 3. Then 
Cre Usres + Crs Ysres + Cro Maas + Crs Mase 
Caran + Cans + Catia + Coetinse = 0 (2). 
Hence, by the antisymmetrical property of Ci, 


Sa ae se (3), 
Ugi29 + thes — tin = 0 
Uyasa — Yarza =0 (4). 


Other instances of case (a) can be obtained by inter- 
changing suffixes that are not already equal, and by turning 
the axes so as to bring 3 into the position of 1, 1 into that 
of 2, and 2 into that of 3. Thus (4) gives 

tyes = Usiee = Uae = Maar = ome = on 5)» 
And also 
thsi = thos = tons = Yin = Ym = Yorn (8), 
Usus = Mane = tons = thar = Mansa = Magn — (7)- 
In case (b), takei = k= m=1,p=2. 
Cy Mane + Creare + Crethare + Crate 
+ Cyathiee + Crotnise + Carta + Continua = 9 (8). 
The last term shows that 
thas = 0 (9), 
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and therefore, by interchange of suffixes, all components of 
class (b) are zero. Also, from the coefficient of cs, 


Use + tae + Mie = 0 (10). 
But in (3) the last term vanishes and we infer 


+ =0 1 
and therefore pues sas 
hare + Usne = 0 12), 

whence, by (10), fe a 
tum = 0 (13). 


‘Thus all components of class (a) are also zero. 
The coefficient of ¢,, in (8) gives 


Un = Une + tha + Yin (14), 


so that the components of class (d) are expressible in terms 
of the three types of class (c). 

No further information is got by transforming com- 
ponents of classes (c) and (¢). Thusifi = k = 1,m= p= 2, 
replacing i or k by j will give a zero component unless j is 
equal to 1; and then the factor ¢,; or ¢,; is zero, and the 
relation holds automatically. Similar considerations apply 
if all of 7, k, m, p are equal. 

We may denote the components of type (5) by A, those 
of type (6) by » + v, and those of type (7) by » — v. Then 
(14) gives 

Urn = Maver = Uases = A+ 2p (15). 
There appear therefore to be three independent: isotropic 
tensors of order 4, obtained by taking each of A, p, v in 
turn equal to 1 and the others to zero. 

In the A tensor, tim» = 1 if i= k and m = p, and in all 
other cases is zero. It is therefore equivalent to 5),5:p, 
which is obviously a tensor of order 4, being the product of 
two tensors of order 2. 

In the p tensor, timp = 1 if i= m, k= p, or if i= p, 
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k=m, andi k. If also i = k, the component is 2. Other 
components are zero. This can be written 

tama = Simin + Bin 3m (16), 
and is obviously a tensor of order 4. 

In the v tensor, tmp = Lifi =m, k = p, and = — lif 
i= -p, k=, and in all other cases is zero. If also i = h, 
Wimp i8 zero. In this case, therefore, 

amy = Bin Sixx — Bp Bim (17). 
This can also be written 
Wamo = €tsn€msp (18), 
for if i= 1, k = 3, erg, = O unless j = 2 and then = 1, But 
then €mjy=1 if m=1, p=3, —1 if m=3, p=1, and 
otherwise = 0, Thus 
thas = 1, tha = — 2 (19), 
with corresponding values for the other components. 


Evidently (17) and (18) represent a tensor of order 4. It 
has already appeared in Chapters 1 and v1. 

The general isotropic tensor of order 4 is therefore 
ABuSmy + H (Bim 3x0 + Bép Spm) + Y (Bem Be9— 94 5em) (20), 
where A, j4, v are scalars. 


EXAMPLE 
Prove that 
Sx 8moin = Smo» 
BimBen + ip 3m) Wie = mp + Mom? 


BimBz9 — 8ip8im) Wie = Pma — Vom 


CHAPTER VIII 
ELASTICITY 


Tn an elastic solid, as in a fluid, the distance between any 
two particles of the body usually varies with the time. 
The body, however, has an equilibrium configuration that 
could persist if the external forces were zero or constant. 
We may take this as a standard of reference. If a 
particle actually at 2, would be at x,— u, in the standard 
configuration, we call u, the displacement at 2,; in practice 
the squares of the u, can usually be neglected. Evidently 
w, is a vector. 

At a point x; + y;, where y;, is small, the displacement is 
u, + v,, where 


ib me bd Ye (1) 


= (Cx — Fu) Ye (2), 


where e,, and &, are the symmetrical and antisymmetrical 
tensors 


cam 4 (5+ 5); ban 4 (We 5) (3). 


If ey, is zero at x,, the displacement has the same form 
as that due to a general displacement u, together with a 
rotation expressed by the tensor &,,. Also, if ey, is every- 
where zero, 


Ba _ 2 (Pmt 
OBtmy iim \O%,— Oatpy 
1 8 (Oty , Ottm 8 (Ory | Atm 
=i Ox, (t+ =) 4 Ox, ( a) 
Berm _ im _ (4). 
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Hence the £, are constants and the rotation is the same 
everywhere. The vanishing of the e, is therefore the neces- 
sary and sufficient condition for a pure rotation. 

Now consider the part of v, due to the em. If we consider 
the quadric surface 
Ca YiYe = 0? (5), 
where 7 is a constant chosen so as to make the surface pass 
through y,, the normal at y; has direction cosines pro- 
portional to x7, and therefore to v;. The displacement 
due to e,, is therefore parallel to the normal at y, to this 
quadric. This coincides with the direction of the radius 
vector to y; if y; is on a principal axis. ‘There are therefore 
three directions such that the relative displacement due 
to the eg is in the direction of y,, and these directions are 
mutually perpendicular. If we take new axes fi fs» & 
along them, the quadric reduces to 

enn'€1? + Cena? + Ca'fs* = 1? (6), 
and all terms é,;' with j # / are zero. The displacement in 
the &, direction is now éy’f,, 80 that all distances in this 
direction are increased in the ratio 1 to 1 + ey’. The dis- 
placement due to the eg, is then the resultant of three 
homogeneous strains parallel to three orthogonal axes. 

We see therefore that the displacement in any small 

neighbourhood can be represented as the combination of 
a rotation with three extensions at right angles. The latter 
express the changes of size and shape of an element of the 
solid. For this reason e,, is called the strain tensor and &_ 
the rotation tensor. Evidently e, has six independent 
components. For we can have 


y= ey, m=%=0 ™, 


making all the e;, zero except @1 = ¢} similarly 2. and e5 
may exist independently of the others. Also, if 


1 = 0, % = ez, = ey (8), 
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all the eg, vanish except é55 imi 
, = =e. Similarl: d 
can be assigned independently. Aaaiiaonh 


In an elastic solid the internal force across an element of 
ae isin general inclined to the surface. If the area of 
e element is dS, the force across it must be specified by 
be components of the form p,,dS parallel to the axes; 
bis n my be regarded as indicating the normal to the 
. If we consider a small parallelepiped with edges 

ou dx,, dx, centred at (2,, %,, 73), the force across the face 
of area dx,dzx, centred at 2, + 4dx, is (py, Py, Drs) A%_ 42%, 
where the py, are evaluated at 2, + 4dz,. The force on the 
opposite face is — (py, Pre, P43) dx,dx, evaluated at 


om +. Of 
2, — }dx,, and the total is Bat daddy. In general the 
force in the x, direction due to the stress across the faces 
of z, constant is Pe de ded, and we take account of 


all faces by using the summation conventi 
; vention. If the 
soenerere of the element has components f, and the 
py orig mass is naeieeitae while if the bodily 
components i 
tions of motion are As a eere ence anes 


phi 2 + 0X, (9). 


The system of quantities py constitutes a symmetri 
tensor of the second order. To prove this, aittae ~e 
plane whose normal has direction cosines a,,, intersecting 
lines through 2, parallel to the coordinate axes at short 
distances from z;; thus a small tetrahedron is formed, 
whose sides are of order J, say. Let the area of the sloping 
face be dS; then those of the others are a,dS. Consider 
now the rate of change of momentum of the matter within 
this tetrahedron. Evidently the contributions from the ac- 
celeration and X, are of the order of the volume, that is, 
of ee The force across dS has magnitude p,,dS. That ities 

yer 
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the face of 2, constant is — p,, times the area of the 
face, that is, — pyid,ydS. But dS is of order J. Hence we 
have 

(Psi — Ges Pes) O (2) = 0 (PY (10), 
and hence if 1 is indefinitely small we have for the stress 
across a plane normal to a,; at 2, 

Ps = Us Pri (11). 

Now consider three perpendicular directions with direction 
cosines ay, (j = 1, 2, 3). The force per unit area across a 
plane perpendicular to one of these axes, in the direction 
of 2, is given by pj. Resolving this along the direction of 
one of the new axes 2, we have, therefore, 


Pir = Ges UrPur = Vs MerPax (12) 
by interchanging i and k; which is precisely the law of 
transformation of a second order tensor. 

Consider again a small parallelepiped centred at x,, with 
edges parallel to the coordinate axes, and form the equa- 
tion of moments about its centre, The contributions from 
fx and X, are of order /* at most, where the edges have 
lengths of order J. The moment about an axis parallel to ~, 
of the stress py in the face x, = constant and parallel to 2, 
is the product of pa into the area of the face and the dis- 
tance of the face from the centre; that is, to order /', 
4pm dar, dx, da. The opposite face makes an equal contribu- 
tion. The stress py, in the face 7, = constant and parallel 
to , tends to turn in the opposite direction. We have, 


therefore, 
(Dar — Piz) 42, datg dey = O (I) (18), 
and therefore when we make J indefinitely small we must 


have 
! Pu = Prs (14), 
and in general 
Pik = Pri (15), 


so that the tensor py, is symmetrical. 
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Now consider the energy interchange between the small 
parallelepiped and its surroundings. The stresses across the 
face of area dx,da, centred at 2, + 4dx, are doing work on 
the element at a rate (py i) dv,dx,, and the contribution 
from the two opposite faces is 


a ' 
an, (Dua tix) dar, daty ders 
Thus in all the stresses are doing work at a rate 
Cc) 
az, (Pestis) dr 


The external forces are doing work at a rate pX,v,dr. The 
kinetic energy of the element is 4pi,*dr, and is increasing 
at a rate pu,f;dr. (We consider the actual specimen of 
matter occupying the element of volume dr at time ¢; thus 
its mass is pdr, and if we keep to the same piece of matter 
at time ¢ + d¢ the mass is unaltered. If we considered the 
variation of energy within a given element of volume we 
should have to allow for the variation of p and the fact that 
the matter moving out of the element is taking its energy 
with it.) The rate of performance of work on the element 
therefore exceeds the rate of increase of kinetic energy by 


a ; 
{2 (Dectis) + pX yy — pid} dr 
a a ea taf 
= {2 (Drstig) + pX pte — ty ( + ex,)} dy, 
by the equations of motion, 
Oui, 
= (puget) dr (16). 

This work is stored up as internal energy of the element of 
volume. Evidently from its form it is a scalar. 


In any elastic solid the internal energy is a definite 
function of the state of the solid. In any change of state 
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the increase of internal energy therefore depends only on 
the initial and final states and not on the method of passage 
from one to the other. Now we have seen that six of the 
are independent, and if the element acquires displace- 
ments 82, in time 5¢ the corresponding increase of internal 
energy is pud (a) dr. There is an apparent asymmetry 
according as i and k are equal or unequal. ‘Thus p,,dr has 
coefficient 8 (2u/x,) = 8e,, but i= 1, k = 2 contributes 
‘pm (2u,/2,), and i= 2, k= 1 contributes p28 (2/2), 
the two together giving 2p.5¢.- But this is the same as 
P23e1 + Pn den, and the whole contribution from the 

of strain is py.5e4d7. Also during the process an 
amount of heat 8Qdr may be absorbed. If then Hdr is the 
internal energy of the element, 

8H = padta + 8Q (17). 

Since B is a definite function of the state of the system, 
and six of the Se, are independent and determine the 
other three, 8 depends on the changes of the temperature 
and of the six independent ¢, and has a definite value in 
whatever order these changes come about. But 


on C0} 
+ a = Pui ( = 2p (18). 


If the absolute temperature is 0, and a certain amount of 
heat 8Q is absorbed without change of any linear dimen- 
sion, the rise of temperature is related to 5Q by the rule 


3Q = pcd0 (19), 


where c is the specific heat at constant strain. If there is 
also a change of strain, since 8Q and 86 are scalars, we must 
have 

8Q = qa5ex + pcdO (20), 


where the q,, constitute a tensor of the second order. But 52 
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and 5Q/0 are perfect differentials. Hence if we replace the 
six independent e,, by e, we can write 


8E = (p, + q,) Bey + pedO (21), 
8Q 1c 
yn be, + G80 (22), 


and 


3 tC) 
Zot a= Zot ani Serta) Ze (60 C3) 


ze:(8)~é0 (6): ao()~a(@) 
It follows at once that if @ is kept constant Zp,Se, and 
Xq,5e, are perfect differentials. Also 


a @ (pc) _ 9 ® (4r\_2ar_ 9 
poet ad= Oa, (%)=e0(4)- 30-6 5) 
and therefore 
ep, 
Gao ag (26). 
If 8@ = 0 we can write 
=p, de, = 8 (27), 


QW = by + 2, + Crelpe, + O(6%) (28). 


The ¢9, Cr, Crs may be functions of @. Then if we retain 
only terms in W up to order e*, 
Pr = Or + Ones (29). 

The c, represent the stresses that would remain if the 
strains e, were removed without change of temperature. 
In most practical cases the original state is one of uniform 
temperature and no stress, so that c, = 0. If there is a rise 
of temperature 6’ under no stress, an element will acquire 
displacements 


where 


0 = On O'Ye (30), 
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where c,, is a second order tensor expressing the thermal 
expansion. Thus 


Cin = ¥ (Cin + ts) Of = Bu! (31), 
where 8, is a symmetrical tensor; and 
Pie = Cte + Cixmp Bp! (32), 


where ¢; my is a fourth order tensor. 

But by hypothesis this deformation takes place under no 
stress and therefore py, = 0. This determines cy, and our 
formula for the stress is 


Pie = Carmo (Cmp — Bmp’) (33). 
The coefficient cym, is the coefficient of eye, in W. 
Since there are six independent e, there are twenty-one 
possible terms in a quadratic form c,,¢,e,, and therefore 
there are twenty-one coefficients ¢;,. They clearly form 
a tensor of order 4; such a tensor in general would have 
eighty-one components, but this satisfies the symmetry 
relations that it is unaltered if we interchange i and k, or 
m and p, or i and & together with m and p together. 
From (26) and (33), 


@ " 
Ya= — Poms FY (Camp) + 02 (Ci, mo Bmp 9") 


== Geno %y (Carmo) + 8CarmsPmp (34) 


if 6’ is small. The second term does not involve the e,,; 
the first is small of the first order in the e,.. 

Many solids are isotropic; that is, they have the same 
properties in all directions. This applies to vitreous (glassy) 
solids and to mixtures of crystals oriented at random. In 
that case a uniform rise of temperature in an element 
gives an equal expansion in all directions and 


v= ab'y, (35) 
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simply; then a is the coefficient of linear expansion and 
Bax On (36). 


‘The second order terms in W, constitute a scalar; and 
we have 


iy 
Ben = Same (37), 


a tensor of order 4, If it is isotropic it must be of the 
form (20) of Chapter vir. Then the linear terms in p,, give 
Pie = Ca,mpemp (38) 
= AB uSmoemo + H (Bim dey + Ber Bem) mp 
+-¥ (8m Sep — Bip Bem) emp 
= ABCmm + H (Bimeme + Sizer) + ¥ (Sim eme— Bsn lxy) 
= ABinemm + H (Cie + Cet) + Y (Cie — Cee) 
= ABlmm + 26 (39), 
the last term vanishing since e,, is symmetrical. 
Then 
2We= Diner (40) 
= NEj¢lmm + 2H OOo 
AAP By (en? + Cea? + Cag? + 2eys? + Leap + 2e49%) (41), 
where A= ey = Ou,/0x, (42). 
The scalars A and » represent properties of the material, 
Both are positive, 
We can also write 
2We= (A+ 2p) A* + Ap (Cg + yy? + Ca" 
— sa — Carr — nem) (43). 
This appears to differ from the form in Love’s Hlasticity, 
1906, p. 100, but the present e, differ from Love’s strain 
components. My ey is the same as his, namely 0u,/22,; but 


my é,, is only half his, so that his assemblage of strain 
components is not a tensor. 
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If all the e,; were equal to one another and therefore to 
3A, we should have 


Wa= (A+ Gp) A* (44), 
Pa=(A+§u)A (i=h); pa=O (b#K) (45). 


In general we write 
A+ Gu=k (46), 


and call & the bulk-modulus. 
2W, = (A+ fp) A? 
+ 2p (1x? + C202 + ya + 2eyq? + Bem? + 2c") 
— Ge (Car® + Can® + C50? + 2ey2€y9 + 2eag@11 + 2611 C22) 
= (A+ fn) A? 
+ 4H (Cu? + e2n* + ean" — C20 a3 — C59 — €11 &2a) 
+ 4p (C2? + a1? + 12) 
= (A+ gu) A* 
+ Gi (C20 — &a)® + (C52 — @11)® + (C11 — @22)® 
+ egg? + Beg? + 669°} (47). 
The coefficient of vanishes if and only if the strain is 
a symmetrical expansion, and may therefore be called the 


distortional strain energy. 
If we allow for variations in temperature, 


Pin = A(A= 800") Bix + 2p Ci» — 0884) 
= {AA — (3A + 22) 28} By + 2nem (48) 


and qa=- 0% is (49). 


Every term in py, is of the first order in the displacements; 
but 20/80 = 1 and therefore gives rise to a constant term. 
This term is 

Jue = 9. Bk dq (50). 
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If8Q=0, so that no heat is lost or gained by conduction, 
pc80 = — qy,5ey, = — 3ka08;,5e, 


= — 3kafdA (51), 
and therefore, if the strain takes place adiabatically, 
CaN (62), 
pe 
as a? 0 


and Pa= (a+ =) Bad + 2ern (53). 
‘Thus in an adiabatic Pe ae the constant A is in- 
creased above its value for a standard disturbance to 0’, 
where 2 

Vedat har (54), 
while jz is unaltered. The bulk-modulus & is therefore also 
increased to k’, where 


= (65). 
Ina simple thermal expansion at zero stress the absorp- 
tion of heat 3Q is equal to pc,50, where c, is called the 
specific heat at zero stress, and is the specific heat mea- 
sured in ordinary experiments. Then 
pty 80 = 8Q = ped0 + qn derx 
= pcd0 + 3ka05y, Se, 
= pc80 + Skal. 3030. 


Bees 
Thus ane (1 + “eg *\ tee (66). 


k 


The equations of motion at constant temperature, if the 
properties \ and y are uniform, can be written 


phi (bad + 2p 0u) + 0X, (67) 


Duy Ot 
“te OA) + Hap, met om) + Xe 


= Le (A+ pw) A}+ nV + PX, (58). 
; 


7 jer 


82 ELASTICITY 


If the changes are adiabatic A must be replaced by 
rv. 

If there is any heat conduction, the absorption of heat 
per unit time by the element of volume dr is 


ze, (K ge) 


where K is the thermal conductivity. Then the equation 
of heat conduction “A 

pe o Pt Wx ba = KV*0 (59). 
If we write 


Ca = 08,0" + Cex! (60), 
so that e,,’ is the strain due to the stresses, 


gust = ok 20° + 300 


a (61), 


and the equation becomes . 
0 28 + Ska X= KV (62). 
The element dr originally had volume 
0 (%— 3 ae Us) ap 
-{i 5 (fat +o met oO (ea) dr (63), 
so that its density was p (1+ A). Hence 


(64). 


This is the equation of continuity. 


CHAPTER IX 
HYDRODYNAMICS 


In comparison with a typical elastic solid, a real fluid shows 
a great resemblance and a fundamental difference. The 
force per unit area across an element of surface parallel to 
a coordinate plane again constitutes a symmetrical tensor 
of order 2, and for the same reasons. The equations of 
motion are still (9) of Chapter vit, and also the rate of 
performance of work on an element of volume dr still has 
the form (16). The difference is that the internal energy in 
a fluid does not depend directly on how much it has been 
deformed. However the fluid is moved about and stirred 
up, provided it returns to its initial position, density, and 
temperature, the initial and final internal energies are 
equal, ‘The deformation, however great it may be, makes 
no contribution; the stresses do work on each element, and 
thereby supply energy, but this is removed in restoring the 
original temperature. If energy of deformation existed in 
a fluid, all particles of it would have a tendency to return 
spontaneously to their original positions when stresses are 
removed, and they have none. Accordingly, while the rate 
of performance of work on the element of volume dr in time 
dt is still p,,0%,/2x,dr, where the symbols have the same 
meanings asinelasticity, wecanno longer assert from this the 
equation (18) of Chapter vim, because the change of internal 
energy is not determinate when the changes of the e,, are 
given. The eg may be as great as we like, but the energy 
does not increase indefinitely apart from changes in density 
and temperature ; and the fluid moves in thesame way under 
the same external forces whatever its previous history. We 
may say that an elastic solid has a memory ; a fluid has none. 
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The stresses in a fluid are related, not to the total de- 
formation, but to the rate of increase of deformation. In 
the former notation these have components 

Oe, Aig , Oty 
“ot 0° 4 (505+ ae) 
The velocities now appear, instead of the displacements 
from an initial configuration, and we now denote the 
velocities, instead of the displacements, by u,. 
We also write 
Buy , Oux\, Quy, — Ou, 

Cn = a(s+ ans fu= a(t - 3) (1), 
so that e,, now denotes the rate of increase of strain and 
£, the local angular velocity, usually called the vorticity. 
Now we say that p;, is linearly related to the new e,,; and 
therefore 

Pix = Pi + Cee, mp&mp (2), 
where P,, is a symmetrical tensor of order 2 that possibly 
does not vanish with the e,,, and ¢¢g, my is a tensor of order 4, 
Further, in a fluid at rest the stress is an isotropic tensor. 
Hence P,, is isotropic and must be of the form — px, 
where p is a scalar, and ¢;p,m, must be of the form (20) 
of Chapter vu. On carrying out the summation with regard 
to m and p the term in v disappears and we have 

Pin = — Pdi + ¥Bix€mm + hese (3), 
where A and p are scalars, at present undetermined. 
Evidently p and A are not both necessary to preserve the 
form, and we may introduce the further convention 

Pau = — 3p (4), 

so that — p is the mean of the three p,, with equal suffixes. 
This gives 

— 8p = — 3p + Bm + Qe (5); 
and therefore A=-—te (6). 

Tf we consider the internal energy as a function of the 
density p and temperature @ alone, we may consider the 
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energy change in a symmetrical expansion. If the density 
increases by §p, there is a contraction in all dimensions in 
the ratio }5p/p and the stresses do work 


~ ta =p? (). 


At the same time there may an absorption or generation 
of heat; then the energy change is 


aE =p? +30 ®) 
and the heat absorbed may be written 
8Q = Mép + pod (9), 
where M is unknown and c is the specific heat at constant 
volume, Then p 
sn = ¢ + i) 3p + pc80 (10), 
8Q MM pe 
g = pt gO qu, 
and the condition that these quantities are perfect differen- 
tials gives A 
p 
m= -05,(2) (12), 
é ae /p 
Fpl) = — 92p(2) (13). 


Tt appears also that p must be a function of p and U) 
alone, for a given material; it does not involve dp/at.* We 
may call p the pressure. 

* ‘This amounts to saying that there is no dissipation of energy in a 
symmetrical compression or expansion. This is true in a gas on the older 
Kinetic theory; but Enskog has shown (Kungl. Svenska Akad. Handlingar, 
63, no. 4, 1922, p. 18) that p can with greater accuracy be given by 

Ro, ndp 
p=9y 9+) a? 
where 7 is a “second” coefficient of viscosity. But n/u is only of the order 
of the square of the ratio of the volume of the molecules themselves to the 
whole volume of the gas. 


In a liquid the coefficient of dp/dt, if it exists, is within the experimental 
error; an analogous statement is true for imperfectly elastic solids. 
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The stress components may now be written 
Pin = — (P+ §plmm) Bix + nee (14), 
where 7 is a function of the density and temperature and 
» expresses a property of the fluid. We call y. the coefficient 


of viscosity. For a uniform fluid we have the equations of 
motion 


afe= mat (P + §H€mm) Bix + » (jul + a} 


+ oe, 
-— o( + Fue) + BVPIs + we ene 
= gg (P— bad) + HVE (15), 
where we now write the scalar 
mm = A (16). 


In time dé the outflow from a volume element dx, da, dx, 
is {0 (pu,)/0x,} drdt. The mass within the element therefore 
decreases at a rate {6 (pu,)/dx,} dr, and we have the equation 
of continuity 

on 
- Zou) (17), 


q, Ou, 
or G7 Pag = 8 (18). 
Here Ly denotes differentiation with regard to the time 


following a given particle of the fluid, so that ae are 
given by 


dx, 
an Us (19), 
d_0,dwya_ 2a C) 
and wi at de dx, Ht Oe, (20). 


In most works on hydrodynamics this operator is denoted 
by D/Dt, but I see no reason for departing from the usual 
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notation for total differentials, since this particular one 
is the only total differential that occurs in hydrodynamics. 
The acceleration components are given similarly by 


Jon Mt un Set (21). 


Consider now the circulation 2 about a closed contour. 
We have 


a5), ty daey = 2 WM ay dec, +f, Uy dey 


-|, (e+ ; ee) dx, + [, udu, (22). 


The integral i u,du, always vanishes ; J X,dz, vanishes 
¢ ec 
for bodily forces derived from a, potential. Then 


ap = | ae (P+ Bnd) Ba + 2p 


mele (p+ Hud) dec + J 7 5h ue) de 


dip + bud) +f 2 ae {u(—26n+ 25) de 


al Ps 
=f paws tua —f2 2 A (ube) de 


2 du a, 2 
aS ad Panel 


If pp = 0 this reduces to the circulation theorem of classical 
hydrodynamics. In many problems of real fluids p is 
small and constant, and A small. Then the integrals 


55 = d (RpA), Joa EWS ins and | % = wd are zero or products 
of two small ped iorte and can be ignored. \3 ~ dp is zero 
cP 


if p is a function of p, which is true in many cases, though 
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exceptions arise when the temperature or the composition 
varies from place to place. In the commonest case, to 
considerable accuracy, 


dQ 
Gi Jame mb de 24), 


and circulation arises only from variations of vorticity in 
the neighbourhood of the contour. It follows that in a 
fluid originally at rest or in irrotational motion circulation 
can arise only through the diffusion of vorticity inwards 
from a boundary. 

The work done on the element dz per unit time exceeds 
the rate of increase of kinetic energy by 


du, 
Pri one dr = Piped (25), 
k 


and Pine = — (P + FHA) Breen + 2eneix 
=— (p+ §pA) A+ 2nenee (26). 


But —pd= 5 ce (27), 


and if p is a function of p oe a constant as a par- 
ticular case) this is a differential with regard to the time 
and yields on the whole zero if the original density is ever 
recovered, This term may therefore be considered to give 
the increase of internal energy due to compression, ‘The 
remainder may be written 


D = — fud?* + 2nenery 
= Gu {eon — esa) ++ (Cas — en) + (Cr — C2)" 
+ Ay (C038 + eay* + €32") (28). 
Thus © is analogous in form to the distortional strain 
energy of elasticity ((47) of Chapter vim). It is essentially 


positive, and therefore represents work done on the fluidand 
continually stored as internal energy. In (8), therefore, if the 
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change is an actual one, the term (27) contributes the p8p/p, 
while (28) contributes to 5Q. So long as the initial and final 
states are given, it is immaterial whether 5Q represents 
heat conducted into the element, absorbed from radiation, 
or generated chemically within it, or mechanical energy 
dissipated into heat by viscosity. 

We notice that ® can vanish if, and only if, 


en = Cra = Ca 
25 = = C12 = 0 (29), 


so that the deformation represents a symmetrical expres- 
sion or contraction. 

If we consider any finite volume, the rate of dissipation 
within it is 


ii @dr= i (— uA? + 2penen) dr (80). 


But eatin = Enbin + out out (31), 
and 
{Ie Sac@n = [ftom eas — fue, (Sn) 


(32), 


where 1, is a direction cosine of the outward normal to the 
boundary. But 


i a fia a 


~ i yw Ends + 4 {I oul ag 


= [fanltnds +a [Eas (39, 


where @/dn denotes differentiation along the outward 
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normal, q is the resultant velocity, and é,, is the vorticity 
vector associated with é,,. Also 


[=f dt r= fh). 


and in all 


[Joe = of 35 as + 4 ff commhénas 


+ 2n {ff (— 48+ faa + ah) ar (88). 


This form is useful in such a problem as that of waves 
on deep water, where the viscosity is small; if it were 
absent we should have a permanent oscillation in a normal 
mode. The vorticity is neaieltle except within a distance 
from the bottom of order (v/y)#, where jp = vp and 2n/y is 
the period of the motion. Then in the second integral £,, 
is zero at the free surface and w, is zero at the bottom, so 
that this integral vanishes, A is negligible everywhere. 
Near the bottom ,, contains a factor proportional to the 
velocity that would exist there on the classical theory, and 
thus can be made indefinitely small for deep water. Hence 
the important term in (35) is 


(fies 
and even this vanishes at the bottom, so that it need only 
be estimated at the free surface. 
If a portion of the fluid is compressed without change of 
temperature, we define the bulk-modulus k by 
oO my (36), 


and if it expands under change of temperature without 
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change of pressure we define a coeflicient of volume 
expansion a by 


14, 
; a =-a (37). 
Then for small changes of temperature and pressure 
sp=p (? = ad?) (38). 
Also in a free expansion under constant pressure 
5Q = Mip + pcs 
a (Pp 
5 {o * (?) op + och 30 (39), 


where the partial differentiation is to be carried out at 
constant density. But this makes 


ih = ak (40), 
and therefore 8Q-= p{e one oh 30 (41), 
so that ¢,, the specific heat at constant pressure, is given by 
ake 
= — 42), 
c+ 73 (42) 


This form is analogous to that of (56) of Chapter vim, the 
present a being the coefficient of volume expansion, and 
the previous one that of linear expansion. 
In an adiabatic change 5Q = 0; then 
— ak *e + pcdd=0 (43), 
and therefore a0Sp = (pc + a®k0) 50 
kA) akd 
c ( & =) = 8s 
so that dp=k ( ae “mye ee (44). 
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Thus the bulk-modulus for adiabatic changes is 


vane (+52) -% (45). 


The equation of heat conduction needs to be modified 
to allow for heat generated internally. In time 5¢, per 
volume dr, we have s B 

8Q = O8t-+ (x 9) at 
= ped + akOASt (46), 


whence pone + akOA = O+ 2(K 22) (47). 


EXAMPLES 


1. Obtain the equations of motion in terms of the stress 
components by considering the momentum of a finite volume of 
any form and applying Green's theorem. 


2. Similarly obtain equation (25). 
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